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Abstract — The global fit of the Standard Model to electroweak precision data, routinely performed by the 
LEP electroweak working group and others, demonstrated impressively the predictive power of electroweak 
unification and quantum loop corrections. We have revisited this fit in view of (i) the development of the 
new generic fitting package, Gfitter, allowing flexible and efficient model testing in high-energy physics, (m) 
the insertion of constraints from direct Higgs searches at LEP and the Tcvatron, and (iii) a more thorough 
statistical interpretation of the results. Gfitter is a modular fitting toolkit, which features predictive the- 
oretical models as independent plugins, and a statistical analysis of the fit results using toy Monte Carlo 
techniques. The state-of-the-art electroweak Standard Model is fully implemented, as well as generic ex- 
tensions to it. Theoretical uncertainties are explicitly included in the fit through scale parameters varying 
within given error ranges. 

This paper introduces the Gfitter project, and presents state-of-the-art results for the global electroweak 
fit in the Standard Model (SM), and for a model with an extended Higgs sector (2HDM). Numerical and 
graphical results for fits with and without including the constraints from the direct Higgs searches at LEP 
and Tevatron are given. Perspectives for future colliders are analysed and discussed. 

In the SM fit including the direct Higgs searches, we find Mh = 116.4 ;j^3^ GeV, and the 2a and 3a allowed 
regions [114, 145] GeV and [[113, 168] and [180, 225]] GeV, respectively. For the strong coupling strength at 
fourth perturbative order we obtain as{M^) = 0.1193 ^q;oq27(6^p) ^ O.OOOl(theo). Finally, for the mass of 
the top quark, excluding the direct measurements, we find mt = 178.2 2 GeV. In the 2HDM we exclude a 
charged-Higgs mass below 240 GeV at 95% confidence level. This limit increases towards larger tan/3, e.g., 
Mh± < 780 GeV is excluded for tan/3 = 70. 
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1 Introduction 

Precision measurements allow us to probe physics at much higher energy scales than the masses 
of the particles directly involved in experimental reactions by exploiting contributions from quan- 
tum loops. These tests do not only require accurate and well understood experimental data but 
also theoretical predictions with controlled uncertainties that match the experimental precision. 
Prominent examples are the LEP precision measurements, which were used in conjunction with the 
Standard Model (SM) to predict via multidimensional parameter fits the mass of the top quark [1], 
prior to its discovery at the Tevatron [2, 3]. Later, when combined with the measured top mass, 
the same approach led to the prediction of a light Higgs boson [4]. Other examples are fits to 
constrain parameters of Supersymmetric or extended Higgs models, using as inputs the anomalous 
magnetic moment of the muon, results on neutral-meson mixing, CP violation, rare loop-induced 
decays of B and K mesons, and the relic matter density of the universe determined from fits of 
cosmological models to data. 

Several theoretical libraries within and beyond the SM have been developed in the past, which, 
tied to a multi-parameter minimisation program, allowed to constrain the unbound parameters of 
the SM [5-8]. However, most of these programs are relatively old, were implemented in outdated 
programming languages, and are difficult to maintain in line with the theoretical and experimental 
progress. It is unsatisfactory to rely on them during the forthcoming era of the Large Hadron 
Collider (LHC) and the preparations for future linear collider projects. Improved measurements 
of important input observables are expected and new observables from discoveries may augment 
the available constraints. None of the previous programs were modular enough to easily allow the 
theoretical predictions to be extended to models beyond the SM, and they are usually tied to a 
particular minimisation package. 

These considerations led to the development of the generic fitting package G fitter [9], designed 
to provide a framework for model testing in high-energy physics. Glitter is implemented in C++ 
and relies on ROOT [10] functionality. Theoretical models are inserted as plugin packages, which 
may be hierarchically organised. Tools for the handling of the data, the fitting, and statistical 
analyses such as toy Monte Carlo sampling are provided by a core package, where theoretical 
errors, correlations, and inter-parameter dependencies are consistently dealt with. The use of 
dynamic parameter caching avoids the recalculation of unchanged results between fit steps, and 
thus significantly reduces the amount of computing time required for a fit. 

The first theoretical framework implemented in Gfitter has been the SM predictions for the elec- 
troweak precision observables measured by the LEP, SLC and the Tevatron experiments. State-of- 
the-art calculations have been used, and - wherever possible - the results have been cross-checked 
against the ZFITTER package [5]. For the W mass and the effective weak mixing angle, which 
exhibit the strongest constraints on the Higgs mass through radiative corrections, the full second 
order corrections are available [11-13]. Furthermore, the corrections of order 0{aa1) and the 
leading three- loop corrections in an expansion of the top-mass-squared (m|) are included. The full 
trhee-loop corrections are known in the large Mh limit, however they turn out to be negligibly 
small [14, 15]. The partial and total widths of the Z are known to leading order, while for the 
second order only the leading corrections are available [16]. Among the new developments in- 
cluded in the SM library is the fourth-order (3NL0) perturbative calculation of the massless QCD 
Adler function [17], contributing to the vector and axial- vector radiator functions in the prediction 
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of the Z hadronic width (and other observables) . It allows to fit the strong coupling constant with 
unique theoretical accuracy [17, 18]. 

Among the experimental precision data used are the Z mass, measured with relative precisions of 
2 • 10~^, the hadronic pole cross section at the Z mass and the leptonic decay width ratio of the 
Z with 10~^ relative precision. The effective weak mixing angle sin^^^g is known from the LEP 
experiments and SLD to a relative precision of 7 • 10~^. The W mass has been measured at LEP 
and the Tevatron to an overall relative precision of 3 • 10~^. The mass of the top quark occurs 
quadratically in loop corrections of many observables. A precision measurement (currently 7-10^^) 
is mandatory. Also required is the precise knowledge of the electromagnetic and weak coupling 
strengths at the appropriate scales. Energy-dependent photon vacuum polarisation contributions 
modify the QED fine structure constant, which at the Z-mass scale has been evaluated to a relative 
precision of 8 • 10~^. The Fermi constant, parametrising the weak coupling strength, is known to 
10~^ relative precision. 

We perform global fits in two versions: the standard ("blue-band") fit makes use of all the available 
information except for the direct Higgs searches performed at LEP and the Tevatron; the complete 
fit uses also the constraints from the direct Higgs searches. Results in this paper are commonly 
derived for both types of fits. 

Several improvements are expected from the LHC [19, 20]. The uncertainty on the VF-boson and 
the top-quark masses should shrink to 1.8 • 10~^ and 5.8 • 10~^ respectively. In addition, the 
Higgs boson should be discovered leaving the SM without an unmeasured parameter (excluding 
here the massive neutrino sector, requiring at least nine additional parameters, which are however 
irrelevant for the results discussed in this paper). The primary focus of the global SM fit would 
then move from parameter estimation to the analysis of the goodness-of-fit with the goal to uncover 
inconsistencies between the model and the data, indicating the presence of new physics. Because 
the Higgs-boson mass enters only logarithmically in the loop corrections, a precision measurement 
is not required for this purpose. Dramatic improvements on SM observables are expected from the 
ILC [21]. The top and Higgs masses may be measured to a relative precision of about 1 • 10~^, 
corresponding to absolute uncertainties of 0.2 GeV and 50 MeV, respectively. Running at lower 
energy with polarised beams, the W mass could be determined to better than 7 • 10~^ relative 
accuracy, and the weak mixing angle to a relative precision of 5 • 10~^. Moreover, new precision 
measurements would enter the fit, namely the two-fermion cross section at higher energies and the 
triple gauge couplings of the electroweak gauge bosons, which are sensitive to models beyond the 
SM. Most importantly, however, both machines are directly sensitive to new phenomena and thus 
either provide additional constraints on fits of new physics models or - if the searches are successful 
- may completely alter our view of the physics at the terascale. The SM will then require extensions, 
the new parameters of which must be determined by a global fit, whose goodness must also be 
probed. To study the impact of the expected experimental improvements on the SM parameter 
determination, we perform fits under the assumption of various prospective setups (LHC, ILC, 
and ILC with GigaZ option). 

As an example for a study beyond the SM we investigate models with an extended Higgs sector of 
two doublets (2HDM). We constrain the mass of the charged Higgs and the ratio of the vacuum 
expectation values of the two Higgs doublets using current measurements of observables from the 
B and K physics sectors and the most recent theoretical 2HDM predictions. 
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The paper is organised as follows. A disquisition of statistical considerations required for the 
interpretation of the fit results is given in Section 2. It is followed in Section 3 by an introduction 
to the Glitter project and toolkit. The calculation of electroweak precision observables, the results 
of the global fit, and its perspectives are described in Section 4. Section 5 discusses results obtained 
for the Two Higgs Doublet Model. Finally, a collection of formulae used in the theoretical libraries 
of Glitter is given in the appendix. We have chosen to give rather exhaustive information here for 
the purpose of clarity and reproducibility of the results presented. 

2 The Statistical Analysis 

The fitting tasks are performed with the Gfitter toolkit described in Section 3. It features the 
minimisation of a test statistics and its interpretation using frequentist statistics. Confidence 
intervals and p- values are obtained with the use of toy Monte Carlo (MC) simulation or probabilistic 
approximations where mandatory due to resource limitations. This section introduces the three 
statistical analyses performed in the paper: (i) determination of SM parameters, (ii) probing the 
overall goodness of the SM, and {in) probing SM extensions and determining its parameters. The 
SM part is represented by the global fit at the electroweak scale (Section 4), while as example for 
beyond SM physics we analyse an extension of the Higgs sector to two scalar doublets (Section 5). 
The statistical treatment of all three analyses relies on a likelihood function formed to measure 
the agreement between data and theory. The statistical discussion below follows in many aspects 
Refs. [22, 23] with additional input from [24, 25] and other statistical literature. 

2.1 Model Parameters 

We consider an analysis involving a set of A'cxp measurements (a;exp)i=i..Arcxp) described by a cor- 
responding set of theoretical expressions (xtheo)i=i..Afcxp- The theoretical expressions are functions 
of a set of A'^mod model parameters (ymod)j=i..Afmod- Their precise definition is irrelevant for the 
present discussion besides the fact that: 

• a subset of (ymod) may be unconstrained parameters of the theory {e.g., the Higgs mass in 
the SM, if the results from the direct searches are not used); 

• another subset of (ymod) are theoretical parameters for which prior knowledge from measure- 
ments or calculations is available and used {e.g., the Z-boson mass and the hadronic vacuum 
polarisation contribution to the running electromagnetic coupling strength); 

• the remaining (ymod) parametrise theoretical uncertainties, which are based on hard-to- 
quantify educated guesswork {e.g., higher order QCD corrections to a truncated perturbative 
series) . 

It may occur that Xcxp or ymod parameters have statistical and theoretical errors, requiring a proper 
treatment for both of these. In the following we use the shorthand notations j/mod {xe:Kp, Xtheo) to 
label both, sets of and individual parameters (measurements, theoretical expressions). 
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2.2 Likelihood Function 



We adopt a least-squares like notation and define the test statistics 




21n£(7/mod) 



(1) 



where the likelihood function, £, is the product of two contributions 




(2) 



The experimental likelihood, >Cexp) measures the agreement between Xtheo and Xexp, while the the- 
oretical likelihood, Z^theo, expresses prior knowledge of some of the Umod parameters. In most cases 
/^exp incorporates well-behaved statistical errors as well as (mostly) non-statistical experimental 
systematic uncertainties. In some instances it may also include theoretical uncertainties and/or 
specific treatments that may account for inconsistent measurements. On the contrary, Ctheo relies 
on educated guesswork, akin to experimental systematic errors, but in most cases less well defined. 
The impact of (mostly strong interactions related) theoretical uncertainties and their treatment 
on the analysis may be strong, as it is the case for the global CKM fit [22, 23]. The statistical 
treatment R&t [22, 23] (described below) is designed to deal with the problem of theoretical er- 
rors in a clear-cut and conservative manner. Evidently though, an ill-defined problem cannot be 
treated rigorously, and results that strongly depend on theory uncertainties must be interpreted 
with care. For the present analysis, by virtue of the large electroweak mass scale so that QCD is in 
the perturbative regime, purely theoretical errors are small and controlled, so that the fit results 
are well behaved. Increasing experimental precision may alter this picture in the future. 

The Experimental Likelihood 

The experimental component of the likelihood is given by the product 



where the iVgxp individual likelihood components Cexpihj) account for observables that may be 
independent or not. The model predictions of the observables depend on a subset of the ymod 
parameters, and are used to constrain those. Ideally, all likelihood components are independent 
(i.e. £exp(i,j) = for i / j) Gaussian functions, each with a standard deviation estimating 
the experimental statistical uncertainty.^ In practise however, one has to deal with correlated 
measurements and with additional experimental and theoretical systematic uncertainties. In ac- 
cordance with the approach adopted by most published analyses, experimental systematic errors 

^ The fitting procedure described in Section 2.3 uses minimisation to obtain the best match between a test 
hypothesis, represented by a certain parameter set, and the data. This requires the use of expected experimentai 
errors corresponding to the test hypothesis in the experimental likehhood, rather than the measured experimental 
errors. However, the expected experimental errors are usually not available for all possible test hypotheses, and the 
measured experimental errors are used instead. This may be a reasonable approximation for test values in close 
vicinity of the measured experimental results. Nonetheless, one should expect that for regions that are strongly 
disfavoured by the likelihood estimator the statistical analysis is less precise, so that large deviations in terms of 
"sigmas" must be interpreted with care. We shall revisit this point in Section 4.2.2 when including results from the 
direct searches for the Higgs boson in the fit. 
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are assumed to express Gaussian standard deviations, so that different systematic errors can be 
added in quadrature.^ Theoretical errors are treated according to the i?fit scheme described below. 

The Theoretical Likelihood 

The theoretical component of the likelihood is given by the product 



The individual components £theo(^) can be constant everywhere in case of no a-priori information, 
be bound, or may express a probabilistic function when such information is reliably available. 
Ideally, one should incorporate in £exp measurements (or equivalent determinations such as Lattice 
gauge theory, provided well-controlled theoretical assumptions are made) from which constraints 
on the ?/mod parameters can be derived. If such constraints are not available, or if a component has 
been explicitly introduced to parametrise theoretical uncertainty, the £theo(0 components must 
be incorporated by hand in Eq. (4). They are statistically ill-defined and can hardly be treated as 
probability density functions. 

In the range fit approach, Rfit, it is proposed that the theoretical likelihoods £theo(0 do not 
contribute to the of the fit when the corresponding Umod parameters take values within allowed 
ranges denoted [j/mod] • Usually these ranges are identified with the intervals [^mod ~ (^theo > ymod + 
o"theo]i where ymod is a best-guess value, and cJtheo is the theoretical systematic error assigned to 
i/mod- Hence all allowed Umod values are treated on equal footing, irrespective of how close they are 
to the edges of the allowed range. Instances where even only one of the 2/mod parameters lies outside 
its nominal range are not considered. This is the unique assumption made in the i?fit scheme: ymod 
parameters for which a-priori information exists are bound to remain within predefined allowed 
ranges. The Rfit scheme departs from a perfect frequentist analysis only because the allowed 
ranges [ymod] do not always extend to the whole physical space. This minimal assumption, is 
nevertheless a strong constraint: all the results obtained should be understood as valid only if all 
the assumed allowed ranges contain the true values of their Umod parameters. Because there is in 
general no guarantee for it being the case, a certain arbitrariness of the results remains and must 
be kept in mind.^ Although in general range errors do not need to be of theoretical origin, but 
could as well parametrise hard-to-assess experimental systematics, or set physical boundaries, we 
will collectively employ the term "theoretical (or theory) errors" to specify range errors throughout 
this paper. 

^This introduces a Bayesian flavour to the statistical analysis. 

^Some j/mod parameters do not have any a-priori information and are hence fuUy unbound in the flt. 

*If a theoretical parameter is bound to an allowed range, and if this range is narrower than what the flt would 
yield as constraint for the parameter if let free to float, the best flt value of this (bound) parameter usually occurs 
on the edge of the allowed range. A modiflcation of this range will thus have immediate consequences for the central 
values of the flt. 
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2.3 Parameter Estimation 

When estimating model parameters one is not interested in the quahty of the agreement between 
data and the theory as a whole. Rather, taking for granted that the theory is correct, one is only 
interested in the quality of the agreement between data and various realisations (models) of the 
theory, specified by distinct sets of 2/mod values. In the following we denote Xmm y ^ absolute 
minimum value of the function of Eq. (1), obtained when letting all the ymod parameters free 
to vary within their respective bounds, with a fit converging at the solution ymod-^ One now 
attempts to estimate confidence intervals for the complete ymod set. This implies the use of the 
offset-corrected test statistics 

Ax^(ymod) = X^(2/mod) - Xmm;y^od ' 

where x^(2/mod) is the for a given set of model parameters ymod- Equation (5) represents the 
logarithm of a profile likelihood. The minimum value Ax^(ymod) is zero, by construction. This 
ensures that, consistent with the assumption that the model is correct, exclusion confidence levels^ 
(CL) equal to zero are obtained when exploring the ymod space. 

In general, the j/mod parameters in Eq. (5) are divided into relevant and irrelevant ones. The 
relevant parameters (denoted a) are scanned for estimation purposes, whereas the irrelevant ones 
(the nuisance parameters //) are adjusted such that Ax^(a, /i) is at a minimum for n = fi. Since in 
frequentist statistics one cannot determine probabilities for certain a values to be true, one must 
derive exclusion CLs. The goal is therefore to set exclusion CLs in the a space irrespective of the 
H values. 

A necessary condition is that the confidence interval (CI) constructed from the Ax^(ymod) test 
statistics provides sufficient coverage, that is, the CI for a parameter under consideration covers 
the true parameter value with a frequency of at least the CL values at the CI boundaries if the 
measurement were repeated many times. For a Gaussian problem, the test statistics follows a x'^ 
distribution [26] and one finds 

1 - CL(a, A) = Prob(Ax^(a, A), dim[a]) , (6) 

where dim[a] is the dimension of the a space, which is the number of degrees of freedom^ of the 
offset-corrected A^^. Here the probability density distribution of A^^ is independent of /i. In a 
non-Gaussian case the CI for a must be evaluated with toy MC simulation for any possible set 
of true /i values using, e.g., a Neyman construction [27] with likelihood-ratio ordering [28, 29].^ 

^The application of the _Rfit scheme in presence of theoretical uncertainties may lead to a non-unique {ymod} 
solution space. 

^Throughout this paper the term confidence level denotes 1 minus the p- value of a given Ax^ (or x^) test statistics, 
and is hence a measure of the exclusion probability of a hypothesis. This is not to be confounded with a confidence 
interval, which expresses an inclusion probability. 

Note that the effective number of degrees of freedom may not always be equal to the dimension of the a space. 
For example, if dim[a] — 2 but a single observable O = /(a) is scanned in a, only one of the two dimensions of o is 
independent, while the other can be derived via O so that the effective dim [a] to be used here is one [23]. Similarly, 
the available observables may only constrain one of the two dimensions of a. Again, the effective dimension to be 
used in Eq. (6) would be one. Intermediate cases, mixing strong and weak constraints in different dimensions of a 
may lead to an ill-posed situation, which can only be resolved by means of a full toy MC analysis. Such an analysis 
is performed at some instances in this paper (see in particular Section 5.2.2 for the two-dimensional case). 

®An ordering scheme is required because the construction of a Neyman CL belt is not unique. It depends on the 
definition of the test statistics used. 
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One may then choose for each a the set of ^ that gives the smallest CL(a). This "supremum" 
approach [24] (also described in Ref. [25] with however a somewhat different meaning) provides 
the most conservative result, which however overcovers in general. (Note also that the approach 
depends on the ordering algorithm used [30]). It may lead to the paradoxical situation that n 
values excluded by the data may be chosen as the true set to determine CL(a). As a modification 
to this scheme, one could only consider fj, values that are within predefined Ax^(a,/i) bounds, 
thus guaranteeing a minimum compatibility with the data [31, 32]. A vast literature on this 
topic is available (see PhyStat conference proceedings and, e.g., Ref. [25]), mostly attempting to 
prescribe a limitation of the space while maintaining good coverage properties.^ We point out 
that the naive "plugin" approach that consists of using the set of fi that minimises Ax^(a,/i) in 
the fit to estimate the true /i is incorrect in general (it is trivially correct if the problem is strictly 
Gaussian, as then the Ax^ distribution is //-independent). It may lead to serious undercoverage 
if the Ax^(a, /x) frequency distribution is strongly dependent on /x (c/. the analysis of the CKM 
phase 7 [24]). 

As a shortcut to avoid the technically challenging full Neyman construction in presence of nuisance 
parameters, one may choose a Gaussian interpretation of the profile likelihood C{a, fi) versus a, 
which corresponds to a MINOS [33] parameter scan. Simple tests suggest satisfying coverage prop- 
erties of the profile likelihood (see, e.g., [34-36]). Mainly because of its simplicity this assumption 
will be adopted for most (though not all) of the results presented in this paper. 

2.4 Probing the Standard Model 

By construction, the parameter estimation via the offset-corrected Ax^ is unable to detect whether 
the SM fails to describe the data. This is because Eq. (5) wipes out the information contained in 
Xmin y d' '^^^^ value is a test statistics for the best possible agreement between data and theory. 
The agreement can be quantified by the p- value 'P(Xmin — Xmin y ^1^^)' which is the tail proba- 
bility to observe a test statistics value as large as or larger than Xmin-gmod ' theory 
underlying the data. It hence quantifies the probability of wrongly rejecting the SM hypothesis. 
In a Gaussian case, y'L;„.,~ can be readily turned into a p-value via Prob(y^- ,nAr,r).^^ In 
presence of non-Gaussian effects, a toy MC simulation must be performed. Again, a full frequentist 
analysis requires the scan of all possible (or "likely") true nuisance parameters, followed by toy 
MC studies to derive the corresponding p-values. Chosen is the set of true ymod that maximises 
'P{Xmin-y J^^)' where here exact coverage is guaranteed by construction (note that in this phase 
no explicit parameter determination is performed so that all ymod are nuisance parameters). 

Such a goodness- of- fit test may not be the most sensitive manner to uncover physics beyond the 
SM (BSM). If the number of degrees of freedom is large in the global fit, and if observables 
that are sensitive to the BSM physics are mixed with insensitive ones, the fluctuations in the 
latter observables dilute the information contained in the global p-value (or deficiencies in the SM 
description may fake presence of new physics). It is therefore mandatory to also probe specific 

'We recall here the reserve expressed in Footnote 1 on page 4 affecting the accuracy of any approach: the 
dependence of the measured errors on the outcome of the observables (determined by a and /i) - if significant - must 
be taken into account. 

^"The corresponding ROOT function is TMath::Prob(...). 
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BSM scenarios. 

2.5 Probing New Physics 

If the above analysis establishes that the SM cannot accommodate the data, that is, the p- value is 
smaller than some critical value, the next step is to probe the BSM physics revealed by the observed 
discrepancy. The goal is akin to the determination of the SM parameters: it is to measure new sets 
of physical parameters y^p that complement the ymod SM parameters. The treatment is identical 
to the one of Section 2.3, using a = {?/np}- Even if the SM cannot be said to be in significant 
disagreement with the data, the estimation of ?/np remains interesting because the most sensitive 
observables, and the precision to be aimed at for their determination can only be derived by this 
type of analysis. Moreover, the specific analysis might be able to faster detect the first signs of a 
discrepancy between data and the SM if the theoretical extension used in the analysis turns out 
to be the right one. 

3 The Gfitter Package 

The generic fitting package Gfitter comprises a statistical framework for model testing and param- 
eter estimation problems. It is specifically designed to provide a modular environment for complex 
fitting tasks, such as the global SM fit to electroweak precision data, and fits beyond the SM. Gfitter 
is also a convenient framework for averaging problems, ranging from simple weighted means using 
or not correlated input data, to more involved problems with non-Gaussian PDFs and/or common 
systematic errors, requiring or not consistent rescaling due to parameter interdependencies. 

Software 

The Gfitter package [37] consists of abstracted object-oriented code in C-I--I-, relying on ROOT 
functionality [10] . The core fitting code and the physics content are organised in separate packages, 
each physics model package can be invoked as a plugin to the framework. The user interfaces Gfitter 
through data cards in XML format, where all the input data and driving options are defined. The 
fits are run alternatively as ROOT macros or executables, interactively or in a batch system. 

Gfitter Parameters and Theories 

Gfitter defines only a single data container, denoted parameter, which can have three distinct 
manifestations according to its use case: 

^^This problem is similar to those occurring in goodness-of-fit (GoF) tests in experimental maximum-likelihood 
analyses. If, for instance, the data sample with respect to which a likelihood analysis is performed is dominated by 
background events with a small but significant signal excess a successful global GoF test would only reveal agreement 
with the background model and say little about the signal. Similarly, a small p-value for the null hypothesis may 
reflect problems in the background description rather than an excess of signal events. A possible remedy here 
would be to restrict the GoF test to signal-like events, or more specifically, to test the GoF in all likelihood bins 
independently. 
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(A) Measurements a^exp that are predicted by the model {e.g., W mass in the SM): parameters 
of this type are not varied in the fit, but contribute to the log-likehhood function through 
comparison between the model prediction and the corresponding measurement. 

(B) Model parameters ymod that are not predicted by the theory but for which a direct mea- 
surement exists {e.g., top mass in the SM): parameters of this type are varied in the fit, and 
they contribute to the log-likelihood function through comparison between the fit parameter 
value and the corresponding measurement. 

(C) Model Parameters ?/mod that are not predicted by the theory and for which no direct mea- 
surement exist {e.g., Higgs mass in the SM), or which parametrise theoretical uncertainties 
according to the -Rfit prescription {of. Section 2.2): parameters of this type are varied freely 
in the fit within bounds (if exist), and they do not contribute themselves to the log-likelihood 
function. 

A parameter is uniquely defined via a name (and optionally an alias to allow the user to declare 
several correlated measurements of the same parameter, and to design theoretical predictions in a 
polymorph class hierarchy) in the data card, and stored in a global parameter container. These 
parameters are objects (of the GParameter class) that cannot be destroyed nor be recreated. Upon 
creation of a parameter, Gfitter searches automatically in the physics libraries for a corresponding 
theory (an object of the GTheory class), identified through the name of the parameter. If a theory 
is found, the corresponding class object is instantiated^^ and the parameter is categorised as of type 
(A); if no theory is found, it is of type (B) or (C) depending on the presence of a measurement in the 
data card.^'^ The categorisation of parameters is performed automatically by Gfitter maintaining 
full transparency for the user. 

Parameter Errors, Ranges, Correlations and Rescaling 

Gfitter distinguishes three types of errors: normal errors following a Gaussian distribution describ- 
ing statistical and experimental systematic errors, a user-defined log-likelihood functions including 
statistical and systematic uncertainties, and allowed ranges describing physical limits or hard-to- 
assess systematic errors (mostly of theoretical origin) . All errors can be asymmetric with respect to 
the central values given. All parameters may have combinations of Gaussian and range errors (but 
only a single user-defined likelihood function). Parameters of type (A) and (B) do not contribute 
to the log-likelihood functions if the theory prediction or fioating parameter value is compatible 
with the central value of the parameter within the ranges of the theoretical errors attributed to 
the parameter (c/. Section 2.1 concerning the implications of the term "theoretical error"). Only 
beyond these ranges, a Gaussian parabolic contribution to the log-likelihood function occurs. For 
example, the combined log- likelihood function of a parameter with central value xq, positive (neg- 
ative) Gaussian error cTQ^^gg ("^Gauss)' ^^'^ positive (negative) theoretical error cr^^^^ ("^theo)' ^ 

^^A GTheory can depend on auxiliary theory objects (derived from GTheory) that are used to outsource complex 
computation tasks. Caching of results from repetitive calculations also benefits from outsourcing. 

^''Measurement results can be given as central value and Gaussian (possibly asymmetric) and/or theoretical errors, 
or as a user-defined log-likelihood function encoded in ROOT objects {e.g. histograms, graphs or functions). 
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given set of ?/mod parameters and theoretical prediction /(ymod) reads^^ 

if: -CTthco ^ /(ymod) -Xo< <Tt"i;eo , 

if: /(ymod) -xo> , 

if: - /(ymod) > C^theo • 

Parameters of type (C) vary freely within the ranges set by the theoretical errors if available, or 
are unbound otherwise. 

Parameters can have correlation coefficients identified and set in the data card via the parameter 
names (and alias if any). These correlations are taken into account in the log- likelihood test 
statistics as well as for the creation of toy MC experiments. 

It is possible to introduce dependencies among parameters, which can be used to parametrise 
correlations due to common systematic errors, or to rescale parameter values and errors with 
newly available results for parameters on which other parameters depend. For example, in the 
global SM fit the experimental value used of the parameter Aa^^^{M'^) depends on UsiM"^). The 
value for as{M'^) used when evaluating Aa^^^{M'^) may have been updated in the meantime, 

or may be updated in each fit step, which leads to a (not necessarily linear) shift of Aa[|^''^(M|) 
and also to a reduced systematic error (for details see Footnote 19 on page 19). The rescaling 
mechanism of Gfitter allows to automatically account for arbitrary functional interdependencies 
between an arbitrary number of parameters. 

Caching 

An important feature of Gfitter is the possibility to cache computation results between fit steps. 
Each parameter holds pointers to the theory objects that depend on it, and the theories keep track 
of all auxiliary theory objects they depend on. Upon computation of the log-likelihood function in 
a new fit step, only those theories (or part of theories) that depend on modified parameters (with 
respect to the previous fit step) are recomputed. More importantly, time intensive calculations 
performed by auxiliary theories that are shared among several theories are made only once per 
fit step. The gain in CPU time of this caching mechanism is substantial, and can reach orders of 
magnitudes in many-parameter fitting problems. 



21og£(y: 



mod ) 



0, 

/(ymod)-(^0+gtteo) 
Cjauss 

/famod)-(3:^0-gt~hco) 



Fitting 

The parameter fitting is transparent with respect to the fitter implementation, which by default 
uses TMinuit [33], but which is extensible via the driving card to the more involved global min- 
ima finders Genetic Algorithm and Simulated Annealing, implemented in the ROOT package 
TMVA [38]. 

^*In the log- likelihood definition of Eq. (7), the central value xo corresponds to the value with the largest likelihood, 
which is not necessarily equal to the arithmetic average in case of asymmetric errors. 
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Parameter Scans and Contours 

Gfitter offers the possibility to study the behaviour of the log-hkelihood test statistics as a function 
of one or two parameters by one- or two-dimensional scans, respectively. If a parameter is of type 
(A), penalty contributions are added to the log- likelihood test statistics forcing the fit to yield the 
parameter value under study. In addition, two-dimensional contour regions of the test statistics 
can be computed using the corresponding TMinuit functionality. 

Toy Monte Carlo Analyses 

Gfitter offers the possibility to perform toy Monte Carlo (MC) analyses repeating the minimisation 
step for input parameter values that are randomly generated around expectation values according 
to specified errors and correlations. For each MC experiment the fit results are recorded allowing 
a statistical analysis, e.g., the determination of a p-value and an overall goodness-of-fit probabil- 
ity. All parameter scans can be optionally performed that way, as opposed to using a Gaussian 
approximation to estimate the p-value for a given scan point (manifestation of true values). 

4 The Standard Model Fit to Electroweak Precision Data 

In recent particle physics history, coined by the success of the electroweak unification and Quantum 
Chromodynamics (QCD), fits to experimental precision data have substantially contributed to our 
knowledge of the Standard Model (SM). The first application of global fits to electroweak data has 
been performed by the LEP Electroweak Working Group [39] in the last decade of the 20*^^ century, 
unifying LEP and SLD precision data. The primary results of these fits were a prediction of the 
top-quark mass (today's fit precision ~ 9 GeV) prior to its discovery, an accurate and theoretically 
well controlled determination of the strong coupling constant at the Z-mass scale (today available 
at the 3NL0 level [17]), and a logarithmic constraint on the Higgs mass establishing that the SM 
Higgs must be light. Other areas related to particle physics where global fits are performed are 
neutrino oscillation [40], leading to constraints on mixing parameters and mass hierarchies, flavour 
physics, with constraints on the parameters of the quark- flavour mixing (CKM) matrix and related 
quantities [22, 41], and cosmology [42], leading to a large number of phenomenological results such 
as the universe's curvature, the relic matter and energy density, neutrino masses and the age of 
the universe. Global fits also exist for models beyond the SM such as Supersymmetry [43, 44] with 
however yet insufficient high-energy data for successfully constraining the parameters of even a 
minimal model so that simplifications are in order. 

We emphasise that the goal of such fits is twofold (c/. Section 2): (i) the determination of the 
free model parameters, and (ii) a goodness-of-fit test measuring the agreement between model and 
data after fit convergence. This latter goal can be only achieved if the model is overconstrained by 
the available measurements. The situation is particularly favourable in the CKM sector, where the 
primary goal of experiments and phenomenological analysis has been moved from CKM parameter 
determination to the detection of new physics via inconsistencies in the CKM phase determination. 
The relatively young field of neutrino oscillation measurements on the contrary does not yet provide 
significant overconstraints of the neutrino flavour mixing matrix. 
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In the following we revisit the global electroweak fit at the Z-mass scale using the Gfitter package. 
We recall the relevant observables, their SM predictions, perform fits under various conditions, 
and discuss the results. 



4.1 Formalism and Observables 

The formal analysis of this section is placed within the framework of the SM. The electroweak 
fit focuses on the parameters directly related to the Z and W boson properties, and to radiative 
corrections to these, providing the sensitivity to heavy particles like the top quark and the Higgs 
boson. The floating parameters of the fit are the Higgs and Z-boson masses, the c, 6, and f-quark 
masses, as well as the electromagnetic and strong coupling strengths at the Z pole. Most of these 
parameters are also directly constrained by measurements included in the fit. 

We have put emphasis on the completeness of the information given in this paper, with a large part 
of the relevant formulae quoted in the main text and the appendices. Readers seeking for a more 
pedagogical introduction are referred to the many excellent reviews on this and related topics (see, 
e.^., Refs. [16, 45-47]). Section 4.1.1 provides a formal introduction of tree-level relations, and 
quantum loop corrections sensitive to particles heavier than the Z . The observables used in the 
global fit and their SM predictions are summarised in Section 4.1.2 and Section 4.1.3 respectively. 
Theoretical uncertainties are discussed in Section 4.1.4. 



4.1.1 Standard Model Tree-Level Relations and Radiative Corrections 

The tree-level vector and axial-vector couplings occurring in the Z boson to fermion-antifermion 
vertex ^/TmIs^j + 5l°/75)/^M are given by^^ 

^(0) _ ^(0) (0) _ ./ ^/ 2 n /ON 

= - 0^°^ - (9) 
9aJ - 9lj 9rj - h ' \^} 

where ff^^Jj) j are the left-handed (right-handed) fermion couplings, and and Jg are respectively 
the charge and the third component of the weak isospin. In the (minimal) SM, containing only 
one Higgs doublet, the weak mixing angle is defined by 

sin2 0^ = l-^. (10) 

z 

Electroweak radiative corrections modify these relations, leading to an effective weak mixing angle 
and effective couplings 



(11) 



9vj = yi4(/|-2Q/sin2^fg) , (12) 

9a 



f = \f7z4^ (13) 



^Throughout this paper the superscript '(0)' is used to label tree-level quantities. 



4.1 Formalism and Observables 



13 



where and are form factors absorbing the radiative corrections. They are given in Eqs. (59) 
and (60) of Appendix A. 3. Due to non-zero absorptive parts in the self-energy and vertex correction 
diagrams, the effective couphngs and the form factors are complex quantities. The observable 
effective mixing angle is given by the real parts of the couplings 

^^l-WsinXV (14) 

Electroweak unification leads to a relation between weak and electromagnetic couplings, which at 
tree level reads 

= „. " ,■■,„..^ . (15) 



^/2(M, 



z 



Radiative corrections are parametrised by multiplying the r.h.s. of Eq. (15) with the form factor 
(1 — Ar)^^. Using Eq. (10) and resolving for Mw gives 



f f 

The form factors p^, and Ar depend nearly quadratically on mt and logarithmically on Mh- 
They have been calculated including two-loop corrections in the on-shell renormalisation scheme 
(OMS) [48-50], except for b quarks where an approximate expression, including the full one-loop 
correction and the known leading two- loop terms oc mf, is provided. The relevant formulae used in 
this analysis are summarised in Appendix A. 3. Since Ar also depends on Mw an iterative method 
is needed to solve Eq. (16). The calculation of Mw has been performed including the complete one- 
loop correction, two- loop and three- loop QCD corrections of order 0{aas) and 0{aa'^), fermionic 
and bosonic two-loop electroweak corrections of order O(a^), and the leading 0{Gl,asm'f) and 
0{G%m^) three-loop contributions [11-13]. Four-loop QCD corrections have been calculated for 
the p-parameter [51-53]. Since they affect the W mass by 2 MeV only, they have been neglected 
in this work. 

For the SM prediction of Mw we use the parametrised formula [11] 

Mw = M'^' - ci dH - C2 dtf + C3 dH^ + C4(dh - 1) - C5 da cq dt 

— C7 dt^ — cs dH dt -|- cg dh dt — cio das + cn dZ , (17) 



with 



lOOGeVy ' VlOOGeVy ' V 174.3 GeV 



91.1875 GeV ' 0.05907 ' 0.119 

where here and below all masses are in units of GeV, and where nit is the top-quark pole mass. Adz 
and Mj{ are the Z and Higgs boson masses, Aa(M^) is the sum of the leptonic and hadronic contri- 
butions to the running QED coupling strength at Af^ (c/. Appendix A.l), as(Af^) is the running 
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strong coupling constant at M| (c/. Appendix A. 2.1), and where the coefficients M^',ci, . . . , cn 
read 



eg 



80.3799 GeV, ci 

0.0000890 GeV, C4 

0.5256 GeV, C7 

0.0000659 GeV, cio 



0.05429 GeV, C2 

0.000161 GeV, C5 

0.0678 GeV, cs 

0.0737 GeV, cn 



0.008939 GeV, 
1.070 GeV, 
0.00179 GeV, 
114.9 GeV. 



The parametrisation reproduces the full result for Miy to better than 0.5 MeV over the range 
10 GeV < Mh < 1 TeV, if all parameters are within their expected (year 2003) 2a intervals [11]. 

The effective weak mixing angle of charged and neutral leptons and light quarks has been com- 
puted [12, 13] with the full electroweak and QCD one-loop and two-loop corrections, and the leading 
three-loop corrections of orders 0{G'^asmf) and ©(G^m^). The corresponding parametrisation 
formula for charged leptons reads 



sin^eis = so + diLn + d2Lj^ + dsL% + d4{Ajj - 1) + d^Aa 

+ deAt + drA'f + dsAtiAn - 1) + dgA^^^ + d^Az , 



(18) 



with 



Lh 
At 



In 



Mh 



100 GeV 



mt 



.178.0 GeV 
and the numerical values 



So 
ds 
de 
dg 



0.2312527, 
3.85 • 10"^ 
-0.002851, 
3.98 • 10""^, 



Aag 



di 
di 
di 
dio 



Mh 



100 GeV 
«s(M|) 
0.117 



4.729 • 10~^ 
-1.85 • 10~^ 
1.82 • 10"^, 
-0.655. 



Az 



1, 



Aa{Mz) 
0.05907 

Mz 
91.1876 GeV 



1, 



d2 = 2.07 • 10"^ 
4 = 0.0207, 
dn = -9.74 • 10" 



Equation (18) reproduces the full expression with maximum (average) deviation of 4.5 • 10~^ 
(1.2 • 10~^), if the Higgs-boson mass lies within 10 GeV < Mh < 1 TeV, and if all parameters are 
within their expected (year 2003) 2a intervals [13]. 

The prediction of the effective weak mixing angle for the remaining light fermions (n, d, s, c quarks 
and neutrinos) differs slightly from the prediction for charged leptons. Again a parametrisation 
formula is provided [13], which is used in this analysis. For bottom quarks, new diagrams with 
additional top-quark propagators enter the calculation and the b quark specific two-loop vertex 
corrections do not exist. Instead we use Eq. (11) and the calculation of (c/. Appendix A. 3), 
which includes the full one- loop correction and the known leading two-loop terms oc mf. 



After completion of this work the two-loop electroweak fermionic corrections to sin^^^g have been published [54]. 
They will be included in future updates of this analysis. 
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4.1 .2 Summary of Electroweak Observables 



The following classes of observables are used in the fit. 
Z resonance parameters: 



Partial Z cross sections: 

Neutral current couplings: 

W boson parameters: 
Higgs boson parameters: 
Additional input parameters: 



Z mass and width, and total e^e Z ^ hadron production 
cross section (i.e., corrected for photon exchange contributions). 

Ratios of leptonic to hadronic, and heavy-flavour hadronic to total 
hadronic cross sections. 

Effective weak mixing angle, and left-right and forward-backward 
asymmetries for universal leptons and heavy quarks. 

W mass and width. 
Higgs mass. 

Heavy-flavour (c, b, t) quark masses (masses of lighter quarks and 
leptons are fixed to their world averages), QED and QCD coupling 
strengths at the Z-mass scale. 



4.1 .3 Theoretical Predictions of Electroweak Observables 



Parity violation in neutral current reactions e~^e~ — > // resulting from the different left and right- 
handed Z-boson couplings to fermions leads to fermion polarisation in the initial and final states 
and thus to observable asymmetry effects. They can be conveniently expressed by the asymmetry 
parameters 



Slj - alj ^ ^ 9v,f/9A,f 



1 + {gvj/gAjf 



(19) 



where only the real parts of the couplings are considered as the asymmetries refer to pure Z 



exchange. For instance, the forward-backward asymmetry = (cr^j — f)/{<^%f + "^s /) 



where the superscript '0' indicates that the observed values have been corrected for radiative 
effects and photon exchange, can be determined from the asymmetry parameters (19) as follows 



l^eAf . 



(20) 



The Af are obtained from Eqs. (19) and (14) using sin^6';fg from the procedure described in the 
previous section. 

Unlike the asymmetry parameters, the partial decay width T j = T{Z ff) is defined inclusively, 
i.e., it contains all real and virtual corrections such that the imaginary parts of the couplings must 
be taken into account. One thus has 



= ANfTo\p{,\{lif 



rUmD + rUmD] , 



(21) 



^^Left-right and forward-backward asymmetries have been also measured for strange quarks, with however insuf- 
ficient precision to be included here. 
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1(3) is the colour factor, Ry{M'^) and ii}j(M|) are radiator functions (defined 



where N[ 
further below), and Fq is given by 







24^/2- 



(22) 



The sin^^^g term entering through the ratio of coupling constants in Eq. (21) is modified by 
the real- valued contribution Ij resulting from the product of two imaginary parts of polarisation 
operators [6] 



where 



«^(Ml)|(l-^Re(4)sin2 0^) . 



The full expression for the partial leptonic width reads [6] 




Ml 




9v,e 



9A,t 



+ 1 



6m^ 
Ml 



1 + 



3a(M| 



vr 



(23) 
(24) 

(25) 



9' 



which includes effects from QED final state radiation. The partial widths for qq final states, T, 
involve radiator functions that describe the final state QED and QCD vector and axial-vector cor- 
rections for quarkonic decay modes. Furthermore, they contain QED QCD and finite quark-mass 
corrections. For the massless perturbative QCD correction, the most recent fourth-order result is 
used [17]. Explicit formulae for the radiator functions are given in Appendix A. 4. The influence 
of non-factorisable EW (8) QCD corrections, Ae^v/qcD) that must be added to the width (21) for 
quark final states is small (less than 10~^). They are assumed to be constant [55, 56], and take 
the values 

—0.113 MeV for u and c quarks, 
^EW/QCD = s —0.160 MeV for d and s quarks, (26) 
-0.040 MeV for the h quark. 

>. 

The total Z width for three light neutrino generations obeys the sum 



Fe + F„ + F, + 3F^ + F 



had 



(27) 



where Fhad = F^ + F^ + Fc + F^ + F;, is the total hadronic Z width. From these the improved 
tree-level total hadronic cross-section at the Z pole is given by 



^0 



127r FgFhad 



Ml 



F2 
^ Z 



(28) 



To reduce systematic uncertainties, the LEP experiments have determined the partial-Z-width 



ratios i?^ 



Fhad/Ff and Ri = F^/Fhadi which are used in the fit. 



The computation of the W boson width is similar to that of the Z boson, but it is only known to 
one electroweak loop. The expression adopted in this analysis can be found in [57]. An improved, 
gauge-independent formulation exists [58], but the difference with respect to the gauge-dependent 
result is small (0.01%) compared to the current experimental error (3%). 
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The value of the QED couphng constant at the Z pole is obtained using three-loop results for 
the leptonic contribution, and the most recent evaluation of the hadronic vacuum polarisation 
contribution for the five quarks lighter than Mz- Perturbative QCD is used for the small top- 
quark contribution. The relevant formulae and references are given in Appendix A.l. 

The evaluation of the running QCD coupling constant uses the known four-loop expansion of 
the QCD /3-function, including three- loop matching at the quark- flavour thresholds (c/. Ap- 
pendix A.2.1). The running of the h and c quark masses is obtained from the corresponding 
four-loop 7-function (c/. Appendix A. 2. 2). All running QCD quantities are evaluated in the 
modified minimal subtraction renormalisation scheme (MS). 

4.1.4 Theoretical Uncertainties 

Within the i?fit scheme, theoretical errors based on educated guesswork are introduced via bound 
theoretical scale parameters in the fit, thus providing a consistent numerical treatment. For exam- 
ple, the effect from a truncated perturbative series is included by adding a deviation parameter, 5th) 
describing the varying perturbative prediction as a function of the contribution from the unknown 
terms. Leaving the deviation parameter floating within estimated ranges allows the fit to adjust 
it when scanning a parameter, such that the likelihood estimator is increased (thus improving the 
fit compatibility). 

f f 

The uncertainties in the form factors and are estimated using different renormalisation 
schemes, and the maximum variations found are assigned as theoretical errors. A detailed nu- 
merical study has been performed in [59] leading to the following real-valued relative theoretical 
errors 

St^py\l-pi\ « 5-10-3^ (29) 
5th4/|l-4l - 5-10-\ (30) 

which vary somewhat depending on the fermion flavour. The corresponding absolute theoretical 
errors are around 2 • 10~^ for both (^th4 ^^'^ ^th^^z treated as fully correlated in the flt. 

These errors, albeit included, have a negligible effect on the fit results. 

More important are theoretical uncertainties affecting directly the M\y and sin'^^gg predictions. 
They arise from three dominant sources of unknown higher-order corrections [11, 13]: (i) 0{a'^as) 
terms beyond the known contribution of 0{Gl,asmf), {ii) 0{a^) electroweak three-loop correc- 
tions, and {Hi) 0(a|) QCD terms. The quadratic sums of the above corrections amount to 

5tY,Mw ~ 4 MeV, (31) 
(^thsin^efg « 4.7•10-^ (32) 

which are the theoretical ranges used in the fit. The empirical W mass parametrisation (17) is 
only valid for a relatively light Higgs boson, Mh ^ 300 GeV, for which the error introduced by 
the approximation is expected to be negligible [11]. For larger Higgs masses, the total theoretical 
error used is linearly increased up to 6i\iM\\r = 6 MeV at Mh = 1 TeV, which is a coarse estimate 
along the theoretical uncertainties given in [11]. 
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Theoretical uncertainties affecting tlie top mass from non-perturbative colour-reconnection effects 
in tlie fragmentation process [60, 61] and due to ambiguities in tlie top-mass definition [62, 63] 
liave been recently estimated to approximately 0.5 GeV each. The systematic error due to shower 
effects may be larger [60]. Especially the colour-reconnection and shower uncertainties, estimated 
by means of a toy model, need to be verified with experimental data and should be included in 
the top-mass result published by the experiments. Both errors have been neglected for the present 
analysis. 

Other theoretical uncertainties are introduced via the evolution of the QED and QCD couplings 
and quark masses, and are discussed in Appendices A.l and A. 2. 

4.2 Global Standard Model Analysis 

The last two decades have been proliferous in providing precision experimental data at the elec- 
troweak scale. Driven by measurements at LEP, SLC and the Tevatron, and significant theoretical 
progress, many phenomenological analyses have been performed, of which we re-examine below the 
global SM fit. The primary goal of this re-analysis is (i) to validate the new fitting toolkit Glitter 
and its SM library with respect to earlier results [5-8], [ii] to include the results from the direct 
Higgs searches at LEP and the Tevatron in the global fit, (Hi) to revisit the impact of theoretical 
uncertainties on the results, and (iv) to perform more complete statistical tests. 

4.2.1 Floating Fit Parameters 

The SM parameters relevant for the global electroweak analysis are the coupling constants of the 
electromagnetic (a), weak (Gp) and strong interactions (as), and the masses of the elementary 
bosons (My, Mz, Mw, Mh) and fermions {mj with / = e, /x, r, z/g, f^, z^r, "w, c, t, d, s, 5,), where 
neutrinos are taken to be massless. The fit simplifies with electroweak unification resulting in a 
massless photon and a relation between the W mass and the electromagnetic coupling a, the Z 
mass, and the weak coupling Gp, according to Eq. (15). Further simplification of the fit arises 
from fixing parameters with insignificant uncertainties compared to the sensitivity of the fit. 

• Compared to Mz the masses of leptons and light quarks are small and/or sufficiently well 
known so that their uncertainties are negligible in the fit. They are fixed to their world 
average values [64]. Only the masses of the heavy quarks, ""^^ rTTc, m;, and rrit, are fioating in 
the fit while being constrained to their experimental values. The top mass uncertainty has 
the strongest impact on the fit. 

• The weak coupling constant Gp has been accurately determined through the measurement 
of the /i lifetime, giving Gp = 1.16637(1) • 10~^ GeV~^ [64]. The parameter is fixed in the 
fit. 

• The leptonic and top-quark vacuum polarisation contributions to the running of the electro- 
magnetic coupling are precisely known or small. Only the hadronic contribution for the five 

^*In the analysis and throughout this paper we use the MS renormahsed masses of the c and b quarks, mc{rnc) 
and mb{rnh), at their proper scales. In the following they are denoted with Trie and mt respectively. 
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lighter quarks, Aa|^J^(M|), adds significant uncertainties and replaces the electromagnetic 
coupling a(M^) as floating parameter in the fit (c/. Appendix A.l). 

With the i?fit treatment of theoretical uncertainties four deviation parameters are introduced in 
the fit. They vary freely within their corresponding error ranges (c/. Section 4.1.4). The theoretical 

uncertainties in the predictions of Mw and sin^^gjj are parametrised by (^th-^w and (5thsin^^gg. 

f f f f 

The form factors and have theoretical errors Sth^z and SthPz: which are treated as fully 

correlated in the fit. 

In summary, the floating parameters in the global electroweak flt are the coupling parameters 
Aa{^^j(M|) and 05 (M^), the masses Mz, rric, f^t and Mh, and four theoretical error param- 
eters. 



4.2.2 Input Data 

A summary of the input data used in the flt is given in the second column of Table 1 , and discussed 
below. 

• The mass and width of the Z boson, the hadronic pole cross section (7^^^, the partial widths 
ratio i?^, and the forward-backward asymmetries for leptons ^fb, have been determined by 
fits to the Z lineshape measured precisely at LEP (see [46] and references therein). Measure- 
ments of the r polarisation at LEP [46] and the left-right asymmetry at SLC [46] have been 
used to determine the lepton asymmetry parameter A(^. The corresponding c and 5-quark 
asymmetry parameters ^c(f))! the forward-backward asymmetries ApB^^\ and the widths ra- 
tios and i?^, have been measured at LEP and SLC [46]. In addition, the forward-backward 
charge asymmetry (Qfb) measurement in inclusive hadronic events at LEP was used to di- 
rectly determine the effective leptonic weak mixing angle sin^^^g [46]. The log- likelihood 
function used in the fit includes the linear correlation coefficients among the Z-lineshape and 
heavy-flavour observables given in Table 2. 

• For the running quark masses rric and mf,, the world average values derived in [64] are used. 
The combined top-quark mass is taken from the Tevatron Electroweak Working Group [65]. 

• For the flve-quark hadronic contribution to a(M^), the most recent phenomenological result 
is used [66] (see also the discussion in [67]). Its dependence on as(M|) requires a proper 
rescaling in the flt (c/. Section 3).^^ 

• The LEP and Tevatron results for the W mass and width are respectively M\y = (80.376 it 
0.033) GeV, Tw = (2. 196 ±0.083) GeV [68], and Mw = (80.432 ±0.039) GeV, Tw = (2.056 ± 
0.062) GeV [69, 70]. Their weighted averages [69], quoted without the correlation coefficient 

^"^ In [66] the light-quark hadronic contribution to a(M|) was found to be Aa[''^'^(i\f|) = 0.02768 ± 0.00022 ± 
0.00002, where the second error singles out the uncertainty from the strong coupling constant for which Qfs(A/|) = 
0.118 ± 0.003 was used. Linear rescaling leads to the modified central value Aa[^^'j(A/|) = 0.02768 + 0.00002 • 
{as{M%)fit — 0.118)/0.003. Since Qs(Mf ) is a free fit parameter and has no uncertainty in a certain fit step the 
error on Aa[^'2j(M|) used in the log-likelihood function does no longer include the contribution from a.s{M%), but 
the corresponding variation is included in the rescaling of the central value only. 



4.2 Global Standard Model Analysis 



20 



Parameter 


Input value 


Free 
in fit 


Results from global EW fits: 
Standard fit Complete fit 


Complete fit w/o 
exp. input in line 


Mz [GeV] 

Tz [GeV] 
<^Ld [nb] 

^FB 

A, w 

Ac 

^6 

.0,c 
^FB 

^FB 

-Rc 

sin2^fff(QFB) 


91.1875 ±0.0021 
2.4952 ± 0.0023 
41.540 ±0.037 
20. 767 ±0.025 
0.0171 ± 0.0010 
0.1499 ±0.0018 

0.670 ±0.027 

qoQ _(_ nof) 

U.C/Z/O _l_ \J.\JjLi\J 

0.0707 ±0.0035 
0.0992 ±0.0016 
0.1721 ±0.0030 
0.21629 ± 0.00066 
0.2324 ±0.0012 
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— 
- 


91.1874 ±0.0021 

2.4959 ±0.0015 
41.477 ±0.014 
20.743 ±0.018 
0.01638 ± 0.0002 

r\ -t A i-r<r^ -1-0 001 1 

0.1478i°:0°JS 

n fiCCO +0.00046 
U.bbS/_Q 0Q045 

0.93470 

0.0741 ±0.0006 
0.1036 ±0.0007 
0.17224 ±0.00006 

n 91 KOI +0.00005 
U.ZIOOI _o. 00007 

0.23143 ±0.00013 


91.1877 ±0.0021 

2.4955 ±0.0015 
41.477 ±0.014 
20.742 ±0.018 
0.01610 ±0.9839 

f-, ^ -1-0 ooos 
n 1 471 tU-uuuo 

-0.0009 

n Rfi»n+o noo32 
U.bb80_o ooo4g 

0.93464 

U.DiOi _o.0005 

inSI +0 0007 
U.lUOl _0.0006 

0.17224 ±0.00006 
0.21580 ±0.00006 

n 9'il 1^1 +0-00012 
U..i01dl _o. 00010 


91.2001 io%°i?^ 
2.4950 ±0.0017 
41.468 ±0.015 
20 71 7 +0 029 

ZU. ( 1 ( _o.025 

0.01616 ±0.0002 

n fiR«n+o noo32 
U.bb80_o 00047 

u.yo'±D'±_o. 00011 

n 07-^7 +0 0004 
U.UIOI _o.0005 

0.1036 ±0.0005 

0.17225 ±0.00006 

0.21580 ±0.00006 
2 3149 +0 00013 

U..i01'iy_o 00009 


Mh [GeV] (°) 


Likelihood ratios 


yes 


on+30[+75] 
°"-23[-41] 


, + 18.3[+28.4] 
110.*_ ^ 3j_ 2 


or,+30[+75] 
°"-23[-41] 


r^y [GeV] 


ou.oyy m u.uzo 
2.098 ±0.048 




80.382 

2.092 tS:""^ 


80.364 ±0.010 
2.091 ±0.001 


an Qi^q+o oio 
su.cioy _o.o2i 

z.uyi _o.oo2 


Trie [GeV] 
mfc [GeV] 

[Hp VI 

a.(M|) 


1.25 ± 0.09 
4.20 ± 0.07 
1 79 4 + 1 9 

2768 ± 22 
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yes 
yes 
yes 


1.25 ±0.09 
4.20 ± 0.07 
172.5 ±1.2 

2772 ± 22 
iiQ2+° °°2® 


1.25 ±0.09 
4.20 ± 0.07 
172.9 ± 1.2 

2767t24 

11 93 +0 0028 


1 7S 9+9-8 

2722^^^ 

1193+00028 

u.iiyo_o 0027 


JthMw [MeV] 
<5thSin20fff (t) 


[-4,4]theo 

[-4.7,4.7]theo 

[^2, 2]thoo 
[^2, 2]thco 
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4 

4.7 

2 
2 


4 
-1.3 
2 
2 





Average of LEP (A« = 0.1465 ± 0.0033) and SLD {At = 0.1513 ± 0.0021) measurements. Tlie complete fit w/o 
the LEP (SLD) measurement gives Ae = 0.14721°;°°°^ {Ae = 0.1463 ± 0.0008 ). (°'ln brackets the 2a errors, ^''''in 
units of 10~®. '^'Rescaled due to as dependency. 



Table 1: Input values and fit results for parameters of the global electroweak fit. The first and second 

columns list respectively the observables/parameters used in the fit, and their experimental values or phe- 
nomenological estimates (see text for references). The subscript "theo" labels theoretical error ranges. The 
third column indicates whether a parameter is floating in the fit. The fourth (fifth) column quotes the 
results of the standard {complete) fit not including (including) the constraints from the direct Higgs searches 
at LEP and Tevatron in the fit. In case of floating parameters the flt results are directly given, while for 
observables, the central values and errors are obtained by individual proflle likelihood scans. The errors are 
derived from the Ax^ profile using a Gaussian approximation. The last column gives the fit results for each 
parameter without using the corresponding experimental constraint in the fit (indirect determination). 
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Mz Vz ol^^ i?° Ali Al^^ Ali A, A, Rl 



Mz 


1 -0.02 


-0.05 


0.03 


0.06 


.0,c 

^FB 


1 0.15 
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-0.02 
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0.07 


Tz 


1 


-0.30 


0.00 


0.00 


.0,6 

^FB 


1 


0.01 


0.06 


0.04 


-0.10 


°had 




1 


0.18 


0.01 


Ac 




1 


0.11 


-0.06 


0.04 








1 


-0.06 


Ab 






1 


0.04 


-0.08 


.0.1 

^FB 








1 










1 


-0.18 



Table 2: Correlation matrices for observables determined by the Z lineshape fit (left), and by heavy 
flavour analyses at the Z pole (right) [46]. 



between mass and width, are used in the fit (c/. Table 1). Since a modest correlation has 
insignificant impact on the fit results^'' it is ignored in the following. 

• The direct searches for the SM Higgs boson at LEP [71] and at the Tevatron [72, 73] use as 
test statistics the negative logarithm of a likelihood ratio, — 21nQ, of the SM Higgs signal 
plus background (s + b) to the background-only (b) hypotheses. This choice guarantees 
— 21n(5 = when there is no experimental sensitivity to a Higgs signal. The corresponding 
one-sided confidence levels CLg+b and CLb describe the probabilities of upward fluctuations 
of the test statistics in presence and absence of a signal (1 — CLb is thus the probability of 
a false discovery). They are derived using toy MC experiments.^^ 

In the modified frequentist approach [74-76], a hypothesis is considered excluded at 95% CL 
if the ratio CLg = CLg+b/CLb is equal or lower than 0.05. The corresponding exclusion 
confidence levels defined by Eq. (6) are given by 1 — CLg and 1 — CLg+b, respectively. The use 
of CLs leads to a more conservative limit [71] than the (usual) approach based on CLg+b-^^ 
Using this method the combination of LEP searches [71] has set the lower limit Mh > 
114.4 GeV at 95% CL. For the Tevatron combination [72, 73], ratios of the 95% CL cross 
section limits to the SM Higgs boson production cross section as a function of the Higgs 
mass are derived, exhibiting a minimum of 1.0 at Mh = 170 GeV. The LEP Higgs Working 
Group provided the observed and expected —2 InQ curves for the s + b and b hypotheses, and 
the corresponding values of the aforementioned confidence levels up to Mh = 120 GeV. The 
Tevatron New Phenomena and Higgs Working Group (TEVNPH) made the same information 
available for 10 discrete data points in the mass range 155 GeV < Mh < 200 GeV based on 
preliminary searches using data samples of up to 3 fb~^ integrated luminosity [73]. For the 
mass range 110 GeV < Mh < 200 GeV, Tevatron results based on 2.4 fb~^ are provided for 

^"A correlation of 0.2 between W mass and width was reported for the Tevatron Run-I resuhs [46]. Assuming the 
same correlation for the LEP and Tevatron combined values of W mass and width leads to an increase of the Xmin 
of the standard fit {complete fit) by 0.09 (0.23). In the complete fit the central value of the Higgs mass estimate is 
unchanged (only the +la bound slightly reduces by 0.6 GeV), whereas a downward shift of 1.1 GeV of the central 
value is observed for the standard fit. In both fits the changes in the other parameters are negligible. 

■^^ For a counting experiment with A'^ observed events and Ns {Ni, ^ Ns) expected signal (background) events, 
one has -InQ = - Nln{Ns/Ni + 1) ~ Ns{l - N/Nt), leading to small -InQ values for large N (signal-like) 
and large — InQ values for small N (background- like). For sufficiently large Ns + Nt, the test statistics — InQ has a 
symmetric Gaussian probability density function. 

■^■^ Assuming a simple counting experiment with a true number of 100 background and 30 signal events, the one- 
sided probability CLg+b to fluctuate to equal or less than 111 observed event is 0.05. The corresponding value 
CLs = 0.05 (which does not represent a probability) is however already reached between 105 and 106 events. 
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— 21nQ [72], however not for the corresponding confidence levels. 

To include the direct Higgs searches in the complete SM fit we interpret the — 21nQ results 
for a given Higgs mass hypothesis^^ as measurements and derive a log-likelihood estimator 
quantifying the deviation of the data from the corresponding SM Higgs expectation. For this 
purpose we transform the one-sided CLg+b into two-sided confidence levels^^ using CLg^^,!,*^'^*^ = 
2CLs+b for CLs+b < 0.5 and CL^;^^'^'^^ ^ 2(1 - CLs+b) for CLs+b > 0.5. The contribution 
to the estimator of the fit is then obtained via = 2 • [Erf^-'^(1 — CL'^I]^^^'^)]'^ , where 
Erf~^ is the inverse error function, and where the underlying probability density function 
has been assumed to be symmetric (c/. Footnote 21 on page 21). 

For the complete mass range available for the LEP searches {Mh < 120 GeV), and for the 
high mass region of the Tevatron searches (155 GeV < Mjj < 200 GeV), we employ the CLg+b 
values determined by the experiments. For the low-mass Tevatron results (110 GeV < Adjj < 
150 GeV), where the CLg+b values are not provided, they are estimated from the measured 
—2lnQ values that are compared with those expected for the s + b hypothesis, and using the 
errors derived by the experiments for the b hypothesis. We have tested this approximation 
in the high-mass region, where the experimental values of CLg+b from the Tevatron are 
provided, and found a systematic overestimation of the contribution to our test statistics 
of about 30%, with small dependence on the Higgs mass. We thus rescale the test statistics 
in the mass region where the CLg+b approximation is used {i.e. 110 GeV < Mh < 150 GeV) 
by the correction factor 0.77.^^ Once made available by the TEVNPH Working Group, this 
approximation will be replaced by the published CLs_(_b values. 

Our method follows the spirit of a global SM fit and takes advantage from downward fluc- 
tuations of the background in the sensitive region to obtain a more restrictive limit on the 
SM Higgs production as is obtained with the modified frequentist approach. The resulting 
curves versus Mh are shown in Fig. 1. The low-mass exclusion is dominated by the LEP 
searches, while the information above 120 GeV is contributed by the Tevatron experiments. 
Following the original figure, the Tevatron measurements have been interpolated by straight 
lines for the purpose of presentation and in the fit which deals with continuous Mh values. 

^■^This procedure only uses the Mh value under consideration, where Higgs-mass hypothesis and measurement are 
compared. It thus neglects that in the SM a given signal hypothesis entails background hypotheses for all Mh valus 
other than the one considered. An analysis accounting for this should provide a statistical comparison of a given 
hypothesis with all available measurements. This however would require to know the correlations among all the 
measurement points (or better: the full experimental likelihood as a function of the Higgs-mass hypothesis), which 
are not provided by the experiments to date. The difference to the hypothesis-only test employed here is expected 
to be small at present, but may become important once an experimental Higgs signal appears, which however has 
insufficient significance yet to claim a discovery (which would allow to discard all other Higgs-mass hypotheses) . We 
thank Bill Murray (RAL) for bringing this point to our attention. 

^''The experiments integrate only the tail towards larger — 21nQ values of the probability density function to 
compute CLs+b (corresponding to a counting experiment with to too few observed events with respect to the s -I- b 
hypothesis), which is later used to derive CLs in the modified frequentist approach. They thus quantify Higgs- 
like (not necessarily SM Higgs) enhancements in the data. In the global SM fit, however, one is interested in the 
compatibility between the SM hypothesis and the experimental data as a whole, and must hence account for any 
deviation, including the tail towards smaller — 21nQ values (corresponding to a counting experiment with too many 
Higgs candidates with respect to the s -I- b hypothesis where, s labels the SM Higgs signal). 

^^The use of Erf ~^ provides a consistent error interpretation when (re)translating the estimator into a confidence 
level via CL = 1 - Prob(x^ 1) = Erf(^P72). 

^^The correction factor reduces the value of the test statistics. As described in Footnote 32, its application has 
little impact on the fit results. 
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Figure 1 : The contribution to the estimator versus Mh derived from the experimental informa- 
tion on direct Higgs boson searches made available by the LEP Higgs Boson and the Tevatron New 
Phenomena and Higgs Boson Working Groups [71-73]. The solid dots indicate the Tevatron mea- 
surements. Following the original figure they have been interpolated by straight lines for the purpose 
of presentation and in the fit. See text for a description of the method applied. 



Constraints on the weak mixing angle can also be derived from atomic parity violation measure- 
ments in caesium, thallium, lead and bismuth. For heavy atoms one determines the weak charge, 
sin 9vir) — A'^. Because the present experimental accuracy of 0.6% (3.2%) for Qw from 
Cs [77, 78] (Tl [79, 80]) is still an order of magnitude away from a competitive constraint on sin^^vi^, 
we do not include it into the fit. (Including it would reduce the error on the fitted Higgs mass 
by 0.2 GeV). Due to the same reason we do not include the parity violation left-right asymmetry 
measurement using fixed target polarised M0ller scattering at low = 0.026 GeV^ [81].^^ 

The NuTeV Collaboration measured ratios of neutral and charged current cross sections in neutrino- 
nucleon scattering at an average ~ 20 GeV^ using both muon neutrino and muon anti-neutrino 
beams [82]. The results derived for the effective weak couplings are not included in this analysis 
because of unclear theoretical uncertainties from QCD effects such as next-to-leading order correc- 
tions and nuclear effects of the bound nucleon parton distribution functions [83] (for reviews see, 
e.g., Refs. [84, 85]). 

Although a large number of precision results for at various scales are available, including recent 
3NL0 determinations at the r-mass scale [17, 18, 86, 87], we do not include these in the fit, because 
- owing to the weak correlation between as(M^) and Mh {cf. Table 3) - the gain in precision 
on the latter quantity is insignificant.^^ Leaving as{M'^) free provides thus an independent and 

^^The main success of this measurement is to have established the running of the weak coupling strength at the 
6.4(7 level. 

^**Including the constraint as(M|) = 0.1212±0.0011 [18] into the fit moves the central value of Mh by -1-0.6 MeV, 
and provides no reduction in the error. 
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theoretically robust determination of the strong coupling at the Z-mass scale. 

The anomaly of the magnetic moment of the muon (g — 2)^ has been measured very accurately to 
a relative precision of 5 • 10^^. Because of the small muon mass the interesting weak corrections 
only set in at a similar size, and this observable is thus not included in the analysis. However, the 
sensitivity of {g — 2)^ to physics beyond the SM (expected to couple to the lepton mass-squared) 
is similar to that of the other observables. 

4.2.3 Fit Results 

All fits discussed in this section minimise the test statistics x^(ymod) defined in Eq. (1). The 
function accounts for the deviations between the observables given in Table 1 and their SM 
predictions (including correlations). Throughout this section we will discuss the results of two fits: 

• The standard ("blue-band") fit, which includes all the observables listed in Table 1, except 
for results from the direct Higgs searches. 

• The complete fit includes also the results from the direct searches for the Higgs boson at LEP 
and the Tevatron using the method described in Section 4.2.2. 

The standard (complete) fit converges at the global minimum value Xmin ~ ^^-^ (Xmin ~ 18.0) for 
13 (14) degrees of freedom, giving the naive p-value Prob(Xmin) 13) = 0.23 (Prob(Xmin' 1^) = 0.21). 
See Section 4.2.5 for a more accurate toy-MC-based determination of the p-value. The results for 
the parameters and observables of the two fits are given in columns four and five of Table 1 together 
with their one standard deviation (a) intervals derived from the Ax^ estimator using a Gaussian 
approximation.^^ We discuss in the following some of the outstanding findings and features of the 
fits. 

Direct and Indirect Determination of Observables, Pulls 

To test the sensitivity of the SM fit to the various input observables, we consecutively disabled 
each of the observables in the fit and performed a log-likelihood scan of the disabled observable. 
The corresponding results and the la intervals are listed in the last column of Table 1. Comparing 
the errors obtained in these indirect determinations with the available measurements reveals their 
importance for the fit. For example, the measurement of Mz is a crucial ingredient, albeit the 
available accuracy is not required. The indirect and direct determinations of Mw are of similar 
precision, such that an improved measurement would immediately impact the fit. The same is 
true for the asymmetry Ai. On the other hand, due to an insufficient precision the heavy quark 
asymmetries Ac and A), do not significantly impact the fit (the fit outperforms the measurements 
by almost two orders of magnitude in precision). 

For further illustration, the pull values obtained from the difference between the fit result and the 
measurement divided by the total experimental error (not including the fit error) are shown for 

^®We have verified tiie Gaussian properties of the fit by sampling toy MC experiments. The results are discussed 
in Section 4.2.4. In the following, unless otherwise stated, confidence levels and error ranges are derived using the 
Gaussian approximation Prob(Ax^, ?idof )• 
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Figure 2: Comparing fit results with direct measurements: pull values for the complete fit (left), and results 
for Mh from the standard fit excluding the respective measurements from the fit (right). 
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Figure 3: Determination of Mh excluding all the sensitive observablcs from the standard fit, except for the 
one given. The results shown are not independent. The information in this figure is complementary to the 
one in the right hand plot of Fig. 2. 
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Table 3: Correlation coefficients between the free fit parameters in the standard fit. The correlations with 
and between the varying theoretical error parameters Sth are negligible in all cases. The correlation between 
Mh and the input parameter Mw amounts to —0.49. 



the complete fit in the left hand plot of Fig. 2 (the standard fit pulls are very similar) . They reflect 
the known tension between the leptonic and hadronic asymmetries, though it is noticeable that no 
single pull value exceeds 3a. The pulls of the c and b quark masses are very small indicating that 
variations of these masses within their respective error estimates has negligible impact on the fit. 
The same observation applies to Mz and Aa[^g^^(M|) (and to a lesser extent even to mt). Thus, 
without significant impact on the goodness-of-fit fit these parameters could have been fixed. 



Correlations 



The correlation coefficients between the fit parameters of the standard fit are given in Table 3. 
Significant are the correlations of —0.40 (+0.31) between In Mh and Aa[j^''^(M^) {mt). An ex- 
cellent precision of these two latter quantities is hence of primary importance for the Higgs-mass 
constraint. The correlation between Aa^^^{M'^) and as(M|) is due to the dependence of the 
hadronic vacuum polarisation contribution on the strong coupling that is known to the fit (c/. 
comment in Footnote 19 on page 19). The correlation coefficients obtained with the complete fit 
are very similar. 



Prediction of thie Higgs Mass 

The primary target of the electroweak fit is the prediction of the Higgs mass. The main results 
are discussed in this paragraph, while more detailed aspects concerning the statistical properties 
of the Higgs mass prediction are presented in Section 4.2.4. The complete fit represents the most 
accurate estimation of Mh considering all available data. We find 

Mh = 116.4 11^3^ GeV (33) 

■^"Fixing rfic, nib, rnt, Mz and /I^ol^I^{Mz) in the fit leads to only an insignificant increase of 0.03 in tlie overall 
Xmini reflecting the little sensitivity of the fit to these parameters varying within the ranges of their (comparably 
small) measurement errors. Of course, this does not prevent Mr to strongly depend on the mt and l\a^^l^{M%) 
input values. 
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where the error accounts for both experimental and theoretical uncertainties. The theory param- 
eters (5th lead to an uncertainty of 8 GeV on Mh, which does however not yet significantly impact 
the error in (33) because of the spread among the input measurements that are sensitive to Mh 
{cf. Fig. 3).^^ As seen in Fig. 12 of Section 4.3, once the measurements are (made) compatible, 
the theoretical errors become visible by the uniform plateau around the Ax^ minimum, and also 
fully contribute to the fit error. The 2a and 3a allowed regions of Mh, including all errors, are 
[114, 145] GeV and [[113, 168] and [180, 225]] GeV, respectively.^^ ^he result for the standard fit 
without the direct Higgs searches is 

Mh = SOtll GeV. (34) 

and the 2a and 3a intervals are respectively [39, 155] GeV and [26, 209] GeV. The 3a upper limit is 
tighter than for the complete fit because of the increase of the best fit value of Mh in the complete 
fit. The contributions from the various measurements to the central value and error of Mh in the 
standard fit are given in the right hand plot of Fig. 2, where all input measurements except for the 
ones listed in a given line are used in the fit. It can be seen that, e.g., the measurements of mt 
and M\Y are essential for an accurate estimation of the Mh- 

Figure 3 gives the complementary information. Among the four observables providing the strongest 
constraint on Mh, namely ^^(LEP), ^^(SLD), ^p'g and Mw, only the one indicated in a given row 
of the plot is included in the fit.^'^ The compatibility among these measurements (c/. Fig. 3) can be 
estimated by (for example) repeating the global fit where the least compatible of the measurements 
(here ^pg) is removed, and by comparing the Xmin estimator obtained in that fit to the one of the 
full fit (here the standard fit). To assign a probability to the observation, the Ax^j^ obtained this 
way must be gauged with toy MC experiments to take into account the "look-elsewhere" effect 
introduced by the explicit selection of the pull outlier. We find that in (1.4 it 0.1)% ("2.5(t") of 
the toy experiments, the Ax^jn found exceeds the Ax^jn = 8.0 observed in current data. 

In spite of the significant anticorrelation between Mh and Aa[j^^^(M^), the present uncertainty 
in the latter quantity does not strongly impact the precision obtained for Mh. Using the theory- 
driven, more precise phenomenological value AajJdl^l) = (277.0± 1.6) • 10"'' [88], we find for the 
standard fit Mh = 80l| GeV. For comparison, with Aa{Jd(^|) = (275.8 ± 3.5) • 10""^ [89], we 
find Mh = 83^26 GeV, reproducing the result form the LEP Electroweak Working Group [39]. 

This is a subtle feature of the -Rfit treatment that we shall illustrate by mean of a simple example. Consider two 
identical uncorrelated measurements of an observable A: 1 ± 1 ± 1 and 1 ± 1 ± 1, where the first errors are statistical 
and the second theoretical. The weighted average of these measurements gives (A) = 1 ± 0.7 ± 1 = 1 ± 1.7, where 
for the last term statistical and theoretical errors (likelihoods) have been combined. If the two measurements only 
barely overlap within their theoretical errors, e.g., 1 ± 1 ± 1 and 3 ± 1 ± 1, their weighted average gives {A) = 2 ± 1. 
Finally, if the two measurements are incompatible, e.g., 1 ± 1 ± 1 and 5 ± 1 ± 1, one finds (A) = 3 ± 0.7, i.e., the 
theoretical errors are only used to increase the global likelihood value of the average, without impacting the error. 
This latter situation occurs in the Mr fits discussed here (although the theoretical errors in these fits are attached 
to the theory predictions rather than to the measurements, which however does not alter the conclusion). 

A fit in which the estimated CLg+b values of the Tevatron searches in the mass region 110 GeV < Mh < 150 GeV 
are not rescaled with the correction factor 0.77 {cf.. Section 4.2.2) leads to a significant increase of the Ax^ value 
only for Mh = 150 GeV. At lower masses the contributions of the direct searches at the Tevatron are small. The 
central value of Mh as well as the la and 3a allowed regions are unchanged; only the 2a intervall is slightly reduced 
to [114, 144] GeV without the correction factor. 

■^■^The uncertainty in the ymod parameters that are correlated to Mh (mainly Aa^^^{Mz) and mt) contributes to 
the errors shown in Fig. 3, and generates a correlations between the four Mh values found. 
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Prediction of thie Top Mass 

Figure 4 shows the Ax^ = ~ Xmin Profile as a function of rrit obtained for the complete fit (soHd 
Une) and the standard fit (dashed hne), both excluding the direct measurement of the top-quark 
mass from the fit. The one, two and three standard deviations from the minimum are indicated 
by the crossings with the corresponding horizontal lines. From the complete fit we find 

mt = 178.21^;^ GeV, (35) 

which, albeit less precise, agrees with the experimental number indicated in Fig. 4 by the dot with 
1(7 error bars (c/. Table 1). The corresponding result for the standard fit is mt = 177.0 t^°o® GeV. 
The insertion of the direct (LEP) Higgs searches leads to a more restrictive constraint towards 
small top-quark masses. Because of the floating Higgs mass, and its positive correlation with mt, 
the profile of the standard fit exhibits an asymmetry (the constraint is less restrictive towards 
larger mt values), which is opposite to the naive expectation from the dominantly quadratic mt 
dependence of the loop corrections. 

Thie Strong and Eiectromagnetic Coupiings 

From the complete fit we find for the strong coupling at the Z-mass scale 

as{Ml) = 0.1193l[J;{5[J^^ ± 0.0001 , (36) 

where the first error is experimental (including also the propagated uncertainties from the errors 
in the c and h quark masses) and the second due to the truncation of the perturbative QCD series. 
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It includes variations of the renormalisation scale between 0.6 Mz < n < 1.3 Mz [18], of massless 
terms of order ag{Mz) and higher, and of quadratic massive terms of order and beyond ag{Mz) 
{cf. Appendix A. 4).^^ Equation (36) represents the theoretically most robust determination of 
as to date. It is in excellent agreement with the recent 3NL0 result from r decays [17, 18], 
as(M|) = 0.1212±0.0005±0.0008±0.0005, where the errors are experimental (first) and theoretical 
(second and third), the latter error being further subdivided into contributions from the prediction 
of the r hadronic width (and spectral moments), and from the evolution to the Z-mass scale. 
Because of their precision, and the almost two orders of magnitude scale difference, the r and 
Z-scale measurements of as represent the best current test of the asymptotic freedom property of 
QCD. 

Finally, the fit result for Aa^^^{M^) without using the constraint from the phenomenological 
analysis in the fit (but including the constraint from the direct Higgs searches, cf. Table 1) precisely 
establishes a running QED coupling, and can be translated into the determination a~^ {Mz )\{it = 
128.99ib0.08. The result is in agreement with the phenomenological value a~^{Mz)\ph = 128.937ib 
0.030 [66]. 



4.2.4 Properties of the Higgs-Mass Constraint 

We proceed with studying the statistical properties of the constraints (33) and (34). Figure 5 
(top) shows the Ax^ profile versus Mh obtained for the standard fit (outermost envelope). Also 
shown is the 95% CL exclusion region obtained from the direct searches at LEP [71]. It exceeds 
the best fit value of the standard fit. The R&t approach provides an inclusive treatment of all 
types of theoretical uncertainties considered in the fit. Fixing the 5th parameters at zero in the fit 
(which is equivalent to ignoring the corresponding theoretical uncertainties) results in a narrower 
log-likelihood curve, with a +0.6 larger global Xmin value, and a shift in Mh at this minimum of 
+2.4 GeV with respect to the result of the standard fit. The difference between the two envelopes 
obtained with freely varying and fixed 5th parameters is highlighted by the shaded band in Fig. 5 
(top). 

In previous electroweak fits [94] theoretical uncertainties were accounted for by independently shifting the 
SM prediction of each affected observable by the size of the estimated theoretical uncertainty, and taking the 
maximum observed cumulative deviation in Mh as theoretical error. The error envelope obtained this way is 
shown in Fig. 6. The dotted curve in the middle of the shaded band is the result of a fit ignoring all theoretical 
uncertainties. The shaded band illustrates the maximum deviations of the Ax^ curves obtained with shifted 
predictions. Including the systematic uncertainties in this way yields a la interval of [55, 122] GeV and 95% 
(99%) CL upper limits of 162 GeV (192 GeV) respectively. For comparison the solid curve in Fig. 6 shows 
the result of the standard fit using the RUt scheme. More detailed studies of systematic theoretical 

^^The uncertainty related to the ambiguity between the use of fixed-order perturbation theory and the so-called 
contour-improved perturbation theory to solve the contour integration of the complex Adler function has been found 
to be very small (3 • 10~^) at the Z-mass scale [18]. 

Another analysis exploiting the r hadronic width and its spectral functions, but using a different set of spectral 
moments than [18], finds as(M|) = 0.1187 + 0.0016 [87]. An analysis of the r hadronic width relying on fixed-order 
perturbation theory finds as{M%) = 0.1180 + 0.0008, where all errors have been added in quadrature [86]. 

■^'^This result is complementary (though more precise) to the LEP measurements of the scale dependence of a 
using, e.g., small and large-angle Bhabha scattering at low energy [90, 91] and high energies [92], respectively, or 
cross section and asymmetry measurements at high energies [93]. 

■^^The inclusion of the theory errors via freely varying parameters (7?fit) leads to a decrease in the global Xmin of 
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Figure 5: Ax^ as a function of Mh for the standard fit (top) and the complete fit (bottom). The solid 
(dashed) hnes give the resuhs when including (ignoring) theoretical errors. The minimum Xmin ^f the fit 
including theoretical errors is used for both curves in each plot to obtain the offset-corrected Ax^- 
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Figure 6: Ax^ versus Mh with an alternative treatment of theory uncertainties [94]. Shown are the 
results of the standard fit ignoring theoretical uncertainties (dotted line), the regions determined from 
the maximum deviation in Ax^ achieved by shifting the SM predictions of all observables according to 1 
"standard deviation" of the various theory uncertainties (shaded band) and for comparison the result of the 
standard fit (solid curve) in which theoretical uncertainties are included in the calculation. 

uncertainties are reported in [59]. 

The Ax^ curve versus Mh for the complete fit is shown in Fig. 5 (bottom). Again the shaded 
band indicates the difference between the two envelopes obtained with freely varying and fixed 
5th parameters, both normalised to the same Xmin (from the fit with free (5th parameters). The 
inclusion of the direct Higgs search results from LEP leads to a strong rise of the Ax^ curve below 
Mh = 115 GeV. The data points from the direct Higgs searches at the Tevatron, available in the 
range 110 GeV < Mh < 200 GeV with linear interpolation between the points, increases the Ax^ 
estimator for Higgs masses above 140 GeV beyond that obtained from the standard fit. 

We have studied the Gaussian (parabolic) properties of the Ax^ estimator to test whether the 
interpretation of the profile likelihood in terms of confidence levels can be simplified. Figure 7 gives 
the 1 — CL derived for Ax^ as a function of the Mh hypothesis for various scenarios: Gaussian 
approximation Prob(Ax^, 1) of the standard fit including theory errors (dashed/red line), Gaussian 
approximation of the standard fit ignoring theory errors, i.e., fixing all 6th parameters at zero 
(solid/black line), and an accurate evaluation using toy MC experiments ignoring theory errors 
(shaded/green area). Also shown is the complete fit result with Gaussian approximation. The toy 
experiments are sampled using as underlying model the best fit parameters (and corresponding 
observables) obtained for each Mh hypothesis. As described in Section 2.4, such a hypothesis is 
incomplete from a frequentist point of view because the true values of the nuisance parameters are 

the fit. Incompatibilities in the input observables (which may be due to statistical fluctuations) thus attenuate the 
numerical efi'ect of the theoretical errors on the fitted parameter (here Mh)- See Footnote 31 on page 27 for an 
illustration of this effect. 
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Figure 7: The 1 — CL function derived from the Ax^ estimator versus the Mh hypothesis (c/. Fig. 5 
(top) for Ax^ versus Mh) for the standard fit. Compared are the Gaussian approximation Prob(Ax^,l) 
for the standard fit with (dashed/red hne) and without theoretical errors (sohd/black line), respectively, to 
an evaluation based on toy MC simulation for which theoretical errors have been ignored. Also given is the 
result using Prob(.) for the complete fit (dotted/blue hne). 

unknown. However, the persuasively Gaussian character of the fit makes us confident that our 
assumption is justified in the present case (c/. the additional discussion and tests in Section 4.2.5). 
The correlations given in Table 2 are taken into account for the generation of the toy experiments. 
Theoretical errors being of non-statistical origin have been excluded from this test, which aims 
at gauging the statistical properties of the test statistics. The curves in Fig. 7 show agreement 
between the Gaussian approximation without theoretical errors, and the toy MC result. It proves 
that the fit is well behaved, and the Ax^ estimator can be interpreted as a true function. 

Figure 8 shows the 68%, 95% and 99% CL contours for the variable pairs rrit vs. Mh (top) 
and Aa^^^{M'^) vs. Mh (bottom), exhibiting the largest correlations in the fits. The contours 
are derived from the Ax^ values found in the profile scans using Prob(Ax^,2) (c/. discussion in 
Section 2.3). Three sets of fits are shown in these plots: the largest /blue (narrower /purple) allowed 
regions are derived from the standard fit excluding (including) the measured values (indicated 
by shaded/light green horizontal bands) for respectively nit and Aa^^^{M'^) in the fits. The 

Examples from other particle physics areas, such as the determination of the CKM phase 7 via direct CP violation 
measurements in B decays involving charm, show that this approximation can lead to severe undercoverage of the 
result [24]. As described in Section 2, the full treatment would require a numerical minimisation of the exclusion 
CLwith respect to any true SM (nuisance) parameter set used to generate the toy MC samples (c/. Refs. [24, 25]). 
More formally, this corresponds to solving CL(AfH) = min^CL^(MH), where jj, are the nuisance parameters of the fit 
and CL^(Afo) = J^'^x'CMff ;data) ^(A;^;2|Mh, ^i)dAx^ and where F(Ax^1Mh,m) is the probability density function 
of Ax^ for true Mh and /i determined from toy MC simulation. 
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Figure 8: Contours of 68%, 95% and 99% CL obtained from scans of fits with fixed variable pairs rrit vs. 
Mh (top) and Aa'^^^^Mz) vs. Mh (bottom). The largest/blue (narrower/purple) ahowed regions are the 
results of the standard fit excluding (including) the measurements of m* (top) and Aa|j^^^(M§) (bottom). 
The narrowest / green areas indicate the constraints obtained for the complete fit including all the available 
data. The horizontal bands indicate the la regions of respectively the mt measurement and Aa^^^^{M'^) 
phenomcnological determination. 
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Figure 9: Contours of 68%, 95% and 99% CL obtained from scans of fits with fixed variable pairs Mw vs. 
mt- The largest /blue allowed regions are the results of the standard fit excluding the measurements of Mw 
and TOj. The narrow/yellow (narrowest /green) areas indicate the constraints obtained for the complete fit 
excluding (including) the corresponding measurements. The horizontal bands indicate the la regions of the 
measurements (world averages). 



correlations seen in these plots are approximately linear for In Mh {cf. Table 3). The third set 
of fits, providing the narrowest constraints, uses the complete fit, i.e., including in addition to 
all available measurements the direct Higgs searches. The structure of allowed areas reflects the 
presence of local minima in the bottom plot of Fig. 5. 

Figure 9 compares the direct measurements of Mw and mt, shown by the shaded/green la 
bands, with the 68%, 95% and 99% CL constraints obtained with again three fit scenarios. The 
largest /blue (narrowest/green) allowed regions are again the result of the standard fit {complete 
fit) excluding (including) the measured values of Mw and mt- The results of the complete fit 
excluding the measured values are illustrated by the narrower /yellow allowed region. The allowed 
regions of the indirect determination is significantly reduced with the insertion of the direct Higgs 
searches. Good agreement is observed between (i) indirect determination without (largest/blue 
area) and with (narrower /yellow area) the direct Higgs searches, and (ii) the direct measurements 
(shaded/green bands). 



4.2.5 Probing the Standard Model 

We evaluate the p- value of the global SM fit following the prescription outlined in Section 2.4. A toy 
MC sample with 10 000 experiments has been generated using as true values for the SM parameters 
the outcomes of the global fit (see the remarks below and in Section 2.4 and Footnote 38 on 
page 32 about the limitation of this method). For each toy simulation, the central values of all the 
observables used in the fit are generated according to Gaussian distributions around their expected 
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Figure 10: Result of the MC toy analysis of the complete fit. Shown are the Xmin distribution of a toy MC 
simulation (open histogram), the corresponding distribution for a complete fit with fixed Sth parameters at 
zero (shaded/green histogram), an ideal distribution assuming a Gaussian case with ridof = 14 (black 
line) and the p- value as a function of the Xmin of the global fit. 



SM values (given the parameter settings) with standard deviations equal to the full experimental 
errors taking into account all correlations.^^ It is assumed that central values and errors are 
independent. The Rfit treatment of theoretical uncertainties allows the fit to adjust theoretical 
predictions and parameters at will within the given error ranges, and - as opposed to measurements 
- the theoretical parameters cannot be described by a probability density distribution and are thus 
not fluctuated in the toy MC. For each toy MC sample, the complete fit is performed (i. e., including 
the results from the direct Higgs searches) yielding the Xmin distribution shown by the light shaded 
histogram in Fig. 10. The distribution obtained when fixing the (^th parameters at zero is shown 
by the dark shaded/green histogram. Including the theoretical uncertainties reduces the number 
of degrees of freedom in the data and hence shifts the distribution to lower values. Overlaid is the 
function expected for Gaussian observables and 14 degrees of freedom. Fair agreement with the 
empirical toy MC distribution for fixed 5th is observed. 

The monotonously decreasing curves in Fig. 10 give the p- value of the SM fit as a function of Xmim 
obtained by integrating the sampled normalised x^ function between Xmin infinity. The value 
of the global SM fit is indicated by the arrow. Including theoretical errors in the fit gives 

p-value(data|SM) = 0.22 ± 0.01 _o.02 , (37) 

where the first error is statistical, determined by the number of toy experiments performed, and 
the second accounts for the shift resulting from fixed parameters. The probability of falsely 

Since only bounds on AIh are available with no probability density information given within these bounds, a 
random generation of Mr toy measurements is not possible. This experimental input is thus kept unchanged for all 
toy MC experiments. 
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Figure 1 1 : P-value of the electroweak fit versus Mh as obtained from toy MC simulation. The error band 
represents the statistical error from the MC sampling. 



rejecting the SM, expressed by the result (37), is sufficient and no significant requirement for 
physics beyond the SM can be inferred from the fit. 

To validate the p^_bcst fit ^ min^p^ assumption used in the above study, we have generated several 
true parameter sets (/i) in the vicinity of the best fit result (varying parameters incoherently by 
±lc7 around their measurement errors), and repeated the toy-MC based p-value evaluation for 
each of them. The probability density distributions derived from these tests have been found 
to be compatible with each other, leading to similar p-values in all cases studied. It supports the 
robustness of the result (37). 

We have extended the above analysis by deriving p-values for the standard fit as a function of the 
true Higgs mass. The results are shown in Fig. 11. For values of Mh around 80 GeV, corresponding 
to the Xmin of standard fit, p-values of about 0.25 are found. '^^ With higher Mh the p-value 
drops reaching the 2(t level at Mh = 190 GeV and the 3(T level at Mh = 270 GeV. 



4.3 Prospects for the LHC and ILC 

The next generation of particle colliders, namely LHC and ILC, have the potential to significantly 
increase the precision of most electroweak observables that are relevant to the fit. This will improve 
the predictive power of the fit, and - in case of a Higgs discovery - its sensitivity to physics beyond 
the SM by directly confronting theory and experiment, and by testing the overall goodness-of-fit 
of the SM. 

''"By fixing Mr the number of degrees of freedom of the fit is increased compared to the standard fit resuhing in 
a larger average Xmin a-nd thus in a larger p-value. 
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At the LHC the masses of the W boson and the top quark are expected to be measured with 
precisions reaching a{Mw) = 15 MeV and a{mt) = 1.0 GeV [19, 20, 95, 96], respectively. At the 
ILC it is expected that the top mass can be measured to an experimental precision of approximately 
cr{mt) = 50 MeV using a threshold scan and an adapted mass definition [21, 99]. This should 
translates into an error of 100-200 MeV on the MS-mass depending on the accuracy of the strong 
coupling constant [21, 99, 100]. More improvements are expected for a linear collider running with 
high luminosity and polarised beams at the Z resonance (GigaZ). The VF-mass can be measured 
to 6 MeV from a scan of the WW threshold [21]. The effective weak mixing angle for leptons 
can be measured to a precision of 1.3 • 10~^ from the left-right asymmetry, Alr [21, 101]. At 
the same time, the ratio of the Z leptonic to hadronic partial decay widths, R^, can be obtained 
to an absolute experimental precision of 0.004 [102]. These numbers do not include theoretical 
uncertainties since it is assumed that substantial theoretical progress will be realised in the years 
left before these measurements are possible. 




Figure 12: Constraints on Mh obtained for the four scenarios given in Table 4, assuming the improvement 
cr(Aa^^'[j(A/^)) — 7 ■ 10^^ for all prospective curves. Shown are, from wider to narrower Ax^ curves: present 
constraint, LHC expectation, ILC expectations with and without GigaZ option. The la errors for Mh given 
in Table 4 correspond to the Ax^ = 1 intervals obtained from these graphs. The shaded bands indicate the 
effect of theoretical uncertainties. 



''^CMS expects a systematic (statistical) precision of better than 20 MeV (10 MeV) for an integrated luminosity 
of 10 fb~^ [20, 97]. It uses a method based on solely the reconstruction of the charged lepton transverse momentum, 
which has reduced systematic uncertainties compared to reconstructing the transverse W mass, with the downside 
of a smaller statistical yield. In an earlier study using the transverse-mass method, ATLAS finds a systematic 
(statistical) uncertainty of better than 25 MeV (2 MeV), for the same integrated luminosity [19]. Combining both, 
lepton channels and experiments, a final uncertainty of about 15 MeV is anticipated in [95], which is used here. A 
recent ATLAS study [98], superseding their previous results, finds that uncertainties of a{Mw) ~ 7 MeV may be 
achievable for each lepton channel (with similar uncertainties for both aforementioned experimental approaches), by 
heavily relying on the calibration of the lepton momenta and reconstruction efficiencies at the Z pole. Using this 
a{Mw) in the fit improves the Mh determination for the LHC prospective from Mh = 120^33 to Mh = I20I26 
(using the improved AaJj^j^(M|) error of 7 ■ 10"^ for both fits, c/. Table 4). 
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Quantity 


Present 


Expected uncertainty 
LHC ILC 


GigaZ (ILC) 


Mw [ MeV] 


25 


15 


15 


6 


mt [ GeV] 


1.2 


1.0 


0.2 


0.1 


sin^efg [10-5] 


17 


17 


17 


1.3 


i?? [10-2] 


2.5 


2.5 


2.5 


0.4 




22 (7) 


22 (7) 


22 (7) 


22 (7) 


Mh{= 120 GeV) [ GeV] 


+ 56 (+52\ r+391 
-40 I- 39; L-31 J 


+ 45 (+A2\ r+301 
-35 l-33j L-25J 


+ 42 r+28i 
-33 L-23J 


+ 27 (+20\ r+81 
-23 l-lSj L-7J 


a,(M|) [10-4] 


28 


28 


27 


6 



Table 4: Measurement prospects at future accelerators for key observables used in the electroweak fit, and 
their impact on the electroweak fit. The columns give, from the left to the right: present errors, the expected 
uncertainties for the LHC with 10 fb^^ integrated luminosity, the ILC without and with the option to run 
at the Z resonance and along the W^-pair production threshold (GigaZ) for one year of nominal running. 
The estimated improvement for Aa[j^j(M^) (given in parenthesis of the corresponding line) over the current 
uncertainty is unrelated to these accelerators, and must come from new low-energy hadronic cross section 
measurements and a more accurate theory (see text). The lower rows give the results obtained for Mh and 

as(M^). For Mh are also given the results with improved Aa|j^'[j(M|) precision (parentheses - this has no 
impact on as{M'^)), and when in addition ignoring the theoretical uncertainties [brackets]. Note that all 
errors obtained on Mh are strongly central value dependent (see text). 

At the time when the new measurements from the LHC experiments, and later the ILC, become 
available, an improved determination of /S.a^^^{M'^) will be needed to fully exploit the new pre- 
cision data. This in turn requires a significant improvement in the quality of the hadronic cross 
section data at energies around the cc resonances and below, and a better knowledge of the c 
and h quark masses entering the perturbative prediction of the cross sections where applicable, 
which serve as input to the dispersion integral. Reference [103] quotes expected uncertainties of 
cj(Aa[f2j(M|)) ~ 7 • 10-5 and 5 • 10-^, compared to presently 22 • 10"^^ if the relative precision 
on the cross sections attains 1% below the J/ip and the T resonances, respectively. The former 
estimate will be used for the present study. Since most of present data is dominated by systematic 
uncertainties, measurements of state-of-the-art experiments with better acceptance and control of 
systematics are needed. High-statistics ISR analyses performed at the B and $ factories already 
provided promising results on many exclusive hadronic channels. New data will also come from 
the BESIII experiment at the BEPCII e'^e~ collider that starts operation in Summer 2008. 

The dominant theoretical uncertainties affecting the electroweak fit arise from the missing higher 
order corrections in the predictions of My^ and sin^^^g (c/. Section 4.1.4), which contribute 
similarly to the error on Mh- They amount to 10 GeV (13 GeV) at Mh = 120 GeV (150 GeV). 
Significant theoretical effort is needed to reduce these. 

A summary of the current and anticipated future uncertainties on the quantities Mw, fnt, sin^^^g, 

i?^, and Aa[f2i(^l)' ^r the LHC, ILC, and the ILC with GigaZ option, is given in Table 4. By 
using these improved measurements the global SM fit (not using the results from direct Higgs 
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searches nor measurements of as{M'^)) results in the constraints on the Higgs mass and as{M'^) 
quoted in Table 4. For all four scenarios the true Higgs mass has been assumed to be Mh = 
120 GeV and the central values for all observables are adjusted such that they are consistent with 
this Mh value. All fits are performed using respectively the present uncertainty on Aa[^^^^(M^), 
and assuming the above-mentioned improvement. For the latter case results for Mh are given 
including (parentheses) and excluding [brackets] theory uncertainties. With the GigaZ option, 
the uncertainty from Aa[^^j(M^) would dominate the overall fit error on Mh if no improvement 
occurred. We emphasise that due (by part) to the logarithmic dependence, the error obtained on 
Mh is strongly Mh dependent: with the same precision on the observables, but central values 
that are consistent with a true value of 150 GeV, one would find Mh = 150 GeV m average, 
i.e., an error increase over the Mh = 120 GeV case of almost 30%. With the GigaZ option and 
the resulting improvement for the uncertainty on as(M^) from the fit is reduced by a factor 
of four. 

The Mh scans obtained for the four scenarios, assuming the improved Aa^^^[M'^) precision to 
be applicable for all future (LHC and beyond) scenarios, are shown in Fig. 12. The shaded bands 
indicate the effect of the current theoretical uncertainties. As expected the theoretical errors 
included with the -Rfit scheme are visible by a broad plateau around the Ax^ minimum. 

A discovery of the Higgs boson at the LHC in the clean decay mode H ^ {H 2i2l') for a 
light (heavy) Higgs would soon allow a precision measurement of Mh beyond the percent level. 
Inserting the measurement into the global electroweak fit would lead to a prediction of the W 
boson mass with 13 MeV error, of which 5 MeV is theoretical. Prediction and measurement could 
be directly confronted. More inclusively, the p- value of the data given the SM could be determined 
as a direct test of the goodness of the SM fit. 

5 Extending the SM Higgs Sector - The Two Higgs Doublet Model 

Two Higgs Doublet Models (2HDM) [104] are simple extensions to the SM Higgs sector, only 
introducing an additional SU{2)l x U{1)y Higgs doublet with hypercharge Y = 1, leading to 
five physical Higgs bosons. Three Higgs bosons {A^, h^, H^) are electrically neutral and the two 
remaining ones (H^) are electrically charged. The free parameters of the 2HDM are the Higgs 
boson masses M^o , Mj^o , Mho and Mh± , the ratio of the vacuum expectation values of the two 
Higgs doublets tan/? = V2/V1, occurring in the mixing of charged and neutral Higgs fields, and the 
angle a, governing the mixing of the neutral CP-even Higgs fields. It should be noted that, in the 
most general 2HDM, tan/3 and hence the corresponding Higgs couplings and mass matrix elements 
depend on the choice of basis for the Higgs fields [105, 106]. 

Models with two Higgs doublets intrinsically fulfil the empirical equality ~ cos^ 6\y . They 
also increase the maximum allowed mass of the lightest neutral Higgs boson for electroweak Baryo- 
genesis scenarios to values not yet excluded by LEP (see, e.g., [107]), and introduce CP violation 
in the Higgs sector. Flavour changing neutral currents can be suppressed with an appropriate 
choice of the Higgs-fermion couplings (see e.g., Ref. [108, 109]). For example, in the Type- 1 2HDM 
this is achieved by letting only one Higgs doublet couple to the fermion sector. In the Type-II 
2HDM [110], which is chosen for this analysis, one Higgs doublet couples to the up-type quarks 
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and leptons only, while the other one couples only to the down-type quarks and leptons. The 
Type-II 2HDM resembles the Higgs sector in the Minimal Supersymmetric Standard Model. It 
fixes the basis of the Higgs fields and promotes tan/J to a physical parameter. 

Our analysis is restricted to the CP conserving 2HDM scalar potential and furthermore we only 
consider observables that are sensitive to corrections from the exchange of a charged Higgs boson. 
In the Type-II 2HDM the charged Higgs-fermion interaction Lagrangian is given by [108, 109] 

Cn^ff = [h+U{MuK (1 - 75) cot/3 + KMd (1 + 75) tan/?)D + h.c.) , (38) 

where U and D are column matrices of three generation up-type and down-type quark fields, 
respectively, Mjj and M£, are the corresponding diagonal mass matrices, and K is the Cabibbo- 
Kobayashi-Maskawa quark-mixing matrix. The charged Higgs interaction has the same structure 
as the charged current mediated by the W . Significant charged Higgs couplings to light quarks 
can occur for large values of tan/J. 

By investigating observables that are sensitive to corrections from a charged Higgs exchange we 
derive constraints on the allowed charged-Higgs mass Mu± and tan/3. Direct searches for the 
charged Higgs have been performed at LEP and the Tevatron. LEP has derived a lower limit of 
Mh± > 78.6 GeV at 95% CL [111], for any value of tan/3. 

5.1 Input Observables 

The constraints on the charged Higgs are currently dominated by indirect measurements, as op- 
posed to direct searches at high-energy accelerators. A multitude of heavy flavour observables 
mainly from S-meson decays is available whose sensitivity to the 2HDM parameters varies how- 
ever substantially, either due to limited experimental precision in case of rare decays, or because 
specific 2HDM contributions are strongly suppressed. The most relevant observables for the search 
of Type-II 2HDM signals are the electroweak precision variable R^, branching fractions of rare 
semileptonic B, D and K decays, and loop- induced radiative B decays. '^^ A summary of the 
experimental input used for this analysis is given in Table 5. 

5.1.1 Hadronic Branching Ratio of Z to b Quarl<s 

The sensitivity of to a charged Higgs boson arises from an exchange diagram modifying the Zbb 
coupling. The corresponding corrections of the SM prediction have been calculated in Ref. [115] 
and are given in Eqs. (6.3) and (6.4) thereof. The left- (right-) handed corrections to the effective 
couplings 5g^^^^ are proportional to cot^/? (tan^/?) and to R/{R — 1) — R\ogR/{R — 1)^, where 
R = m'l/M'jjj^ . The charged-Higgs exchange leads to a decrease of R'^. Neutral Higgs contributions 
can be neglected for small tan/3. For the SM prediction we use the result from the complete 
electroweak fit, -RfJgM = 0.21580 ±0.00006, where the direct measurement of R'^ has been excluded 

''■^Decays of r and fj, leptons can also occur through charged-Higgs tree diagrams giving anomalous contributions 
to the decay parameters {Michel parameters [112]) measured in these decays. Their present sensitivity is however 
not competitive with the other observables (a 95% CL limit of Mjf± > 1.9 GeV • tan/3 is currently achieved from r 
decays [113], see also [114] for a review of the /i decay parameters). 
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(c/. Table 1). It is confronted in the fit with the experimental value R^^^^ = 0.21629 it 0.00066, 
obtained at LEP [46], giving Ai?° = R^^^p - R^^^ = 0.00049 ± 0.00066.' 

5.1.2 The Decay B — ^ Xs7 

The decay B — > X^^ is an effective flavour changing neutral current process occurring only at loop- 
level in the SM. The SM prediction for its branching fraction (B) at NNLO accuracy is (3.15 it 
0.23) • 10"'^ [116], where the theoretical uncertainty is estimated by studying (in decreasing order of 
importance) nonperturbative, parametric, higher-order and rric interpolation ambiguity effects, and 
where all errors have been added in quadrature. Averaging branching fraction measurements from 
the BABAR, Belle and CLEO Collaborations gives B{B Xs-f) = (3.52 ±0.23 ±0.09) • 10"^ [117], 
where the first error is experimental and the second stems from the modelling of the photon energy 
spectrum. The improved NNLO calculation yields a branching fraction approximately 1.5 a lower 
than the NLO calculation [118], resulting in a small tension with the experimental average and 
thus leading to less stringent constraints on the charged Higgs mass. The 2HDM contribution to 
the B{B — > Xgj) arises from a charged Higgs replacing the in the loop from which the photon 
is radiated and is always positive in the type-II model. For the prediction of B{B Xsj) in 
the 2HDM we have used parametrised formulae [119] reproducing the result of [116] within 0.2%. 
While the value of the branching fraction changes with Mjj± and, to a lesser extent with tan/3, the 
associated theoretical uncertainty stays to good approximation constant at 7%. Since it has been 
derived by quadratically combining several error estimates, we treat it as an additional Gaussian 
systematic error in the fit. 

5.1 .3 Leptonic Decays of Charged Pseudoscalar Mesons 

In the SM the leptonic decay of charged pseudoscalar mesons proceeds via the annihilation of the 
heavy meson into a W boson and its subsequent leptonic decay. Angular momentum conservation 
leads to a helicity suppression factor that is squared in the lepton mass. Competitive contributions 
from the charged Higgs sector can therefore occur. Neglecting photon radiation, the leptonic decay 
rate of a pseudoscalar meson P has the form 

r (P ^ £u) = = —fpm^mp 1 2" 1^9192! > (39) 

Tp Sir \ nip J 

where mp {mi) is the mass of the pseudoscalar meson (lepton), jV^igjI the magnitude of the 
CKM matrix element of the constituent quarks in P, and fp is the weak decay constant. 

For P = 5 (implying P± = B^) we use [117] tb± = (1.639±0.009)ps and |T4fe| = (3.81±0.47)-10-3, 
where the latter result has been averaged over inclusive and exclusive measurements. For the B 
decay constant we use the value fp = (216 ± 22) MeV, obtained by the HPQCD Collaboration 
from unquenched Lattice QCD calculations [120]. For meson and lepton masses we use the values 
of Ref. [64]. With these inputs, we find the SM predictions B{B tv) = 1.53 Ig'J^ • 10"^ and 

B{B fiu) = omtoll • 10"^. 

An alternative approach uses for the r.h.s. of Eq. (39) additional constraints from the global 
CKM fit enhancing the information on \Vub\ beyond that of the direct measurement through the 
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fit of the Wolfenstein parameters p, r], and on tlirougli the measurement of the B^B^ mixing 
frequency. This assumes that the measurements entering the fit are free from significant new 
physics contributions. It is certainly the case for the charged Higgs, but cannot be excluded for 
the CP-violation and neutral-i? mixing observables. Hence, albeit using the global CKM fit is 
an interesting test, it cannot replace the direct SM prediction of Eq. (39) based on tree-level 
quantities and lattice calculations only. Not using the direct measurements, the global CKM fit 
gives \Vub\ = (3.44l°-j2) . IQ-^, and for the complete prediction B{B tv) = 0.83 [121]. 
This latter result is about l.Ocr below the one from the "tree-level" determination, and a similar 
discrepancy is found for B ^ pv (cf. Table 5). 

The charged-Higgs amplitude contributes to the leptonic decays modifying Eq. (39) by a scaling 
factor rn- In the Type-II 2HDM the b quark couples only to one of the Higgs doublets at tree 
level so that the scaling factor for the decays B ^ tv and B ^ fiu reads [122] 



1 — mi 



tan/? \ 



(40) 



which can lead to both, an increase and a decrease in the branching fraction, depending on whether 
the and amplitudes interfere constructively or destructively. 

The rare leptonic decay B ^ tu has been observed by the BABAR and Belle Collaborations [123- 
125], with an average branching fraction^^ ^ ^ 5;^ _^ o.33) • 10"''. Only upper limits 

are available for the muon channel so far, the tightest one, B{B fiu) < 1.3 • 10"® at 90% CL 
(—12 lb 20 fitted events), being found by BABAR [127]. For lack of an experimental likelihood we 
use the measured branching fraction of (—5.7 i 7.1stat i 6.8syst 

) X 10"^ 

For P = K, contributions from a charged Higgs are suppressed by {rriK /^b)'^ relative to leptonic 
B decays. Moreover, due to the smaller phase space for hadronic final states, leptonic decays have 
large branching fractions, which - on the other hand - have been measured to an excellent 0.2% 
relative accuracy for i = fi. We follow the approach of Ref. [128] and compare \Vus\ determined 
from helicity suppressed K ^ decays and helicity allowed K njjLU decays, considering the 
expression 

Vus{K ^ ^jiu) Krf(0+ ^ 0+ 



-^£23 



(41) 



which in the SM is equal to 1. The ratio 13{K ^u)/B{Ti fiu) ~ (yusfK)/iyudfn) is used to 
reduce the theoretical uncertainties from the kaon decay constant fx, and from electromagnetic 
corrections in the decay K iiv [128]. The dominant uncertainty in Vus from K — > tt/xi^ decays 
stems from the — > vr vector form factor at zero momentum transfer, /+(0), while V^d determined 
from super-allowed nuclear beta-decays (0^ 0^) is known with very high precision [129]. 

In the 2HDM of Type-II the dependence of -Rf23 due to charged Higgs exchange is given by [128] 



R 



H 



m 



H+ 



(42) 



Updated results from BABAR and Belle have been presented at the recent workshops CKM 2008 and Tau 
2008 [126], leading to the new average B{B Tu) = (1.73 ± 0.35) • 10"". They will be included in future updates 
of this analysis. 



5.2 Results and Discussion 



43 



Parameter 


Experimental value 


Ref. 


SM prediction 


Ref. 


Rb 


U.ziDZy ± U.UUUDO 


[4bJ 


U.ZiDoU it U.UUUUO 


This work 


[5{B XsJ) [10 \ 


o.bz ± u.zo ± u.uy 


[117J 


o.lo ± U.ZO 


[116] 




1.51 ±0.33 


[125] 


1 53"^nQQ ffe IKhll 

-*^*'-'^ — 0.38 L./-^' 1 ^^^iolJ 
0.83 +2;^^ [CKM fit] 


This work 
[121] 


B{B /ii/) [10-^] 


-5.7 ±6.8 ±7.1 


[127] 


0.69l°-ji [/s,|K,|] 
0-37t2;o4 [CKM fit] 


This work 
[121] 




0.42 ±0.12 ±0.05 


[131] 


0.28 ±0.02 


[130] 




1.004 ±0.007 


[128] 


1 





Table 5: Experimental results and SM predictions for the input observables used in the analysis of the 
charged-Higgs sector of the Type-II 2HDM. 



where we use md/ms = 19.5 ± 2.5 [64]. Experimentally, a value of R^23 — 1-004 ± 0.007 is 
found [128], where {fK/fn)/f+{0) has been taken from lattice calculations. It dominates the 
uncertainty on Re23- 



5.1.4 The Semileptonic Decay B Dtu 

Similar to the B ^ decay, the semileptonic decay B Dtv can be mediated a by charged 
Higgs. We follow the arguments of Ref. [130] and use the ratio Rdt/e = B{B — > Dtv) / B{B Dev) 
to reduce theoretical uncertainties from hadronic form factors occurring in the predictions of the 
individual branching fractions. In the Type-II 2HD1VI the ratio Rdt/s can be expressed in the 
following compact form [130] 

R^^^^ = (0.28 ± 0.02) • [l + (1.38 ± 0.03) • Re(C5p) ± (0.88 ± 0.02) • jC^pl^l , (43) 

where C^p = — mfom,- tan^/3/m^± . As for leptonic decays the 2HDM contribution can either lead 
to an increase or decrease in the branching fraction. Equation 43 is the result of an integration 
of the partial width dT{B — > Dii')/dw, assuming no Higgs contribution to — > Deu, and where 
w = vbvd with vb (vd) being the four- velocity of the B (D) meson. 

The ratio of branching fractions has been measured by BABAR to be R^jy^j^ = 0.42 ± 0.12stat ± 
0.05,yst [131]. 



5.2 Results and Discussion 

The theoretical predictions of the Type-II 2IIDIVI for the various observables sensitive to corrections 
from the exchange of charged Higgs bosons have been implemented in a separate library integrated 
as a plug-in into the Glitter framework. Exclusion confidence levels have been derived in two ways: 
(i) for each observable separately, and {ii) in a combined fit. 
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Figure 13: Two-sided 68%, 95% and 99% CL exclusion regions obtained for the various observables (see 
text) in the 2HDM parameter plane M^± versus tan/3. 
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5.2.1 Separate Constraints from Individual Observables 

Constraints in the two-dimensional model parameter plane (tanP, Mfj±) have been derived using 
the individual experimental measurements and the corresponding theoretical predictions of the 
Type-II 2HDM. Figure 13 displays the resulting two-sided 68% (yellow/light), 95% (orange) and 
99% CL (red/dark) excluded regions separately for each of the observables given in Table 5. The 
confidence levels are derived assuming Gaussian behaviour of the test statistics, and using one 
degree of freedom (c/. discussion in Footnote 7 on page 6) , i.e., Prob(Ax^,l). Also indicated in 
the plots is the 95% CL exclusion limit resulting from the direct searches for a charged Higgs at 
LEP [111] (hatched area). 

The figures show that Rf, is mainly sensitive to tan/3 excluding small values (below ~1). B{B — > 
Xg'y) is only sensitive to tan/3 for values below ~1. For larger tan/3 it provides an almost constant 
area of exclusion of a charged Higgs lighter than ~260 GeV. (All exclusions at 95% CL). The 
leptonic observables lead to triangle-shaped excluded areas in the region of large tan/3 and small 
mf{± values. 13{B — > tu) gives the strongest constraint. '^^ For these observables the 2HDM 
contribution can be either positive or negative, because magnitudes of signed terms occur in the 
predictions of the branching fractions giving a two-fold ambiguity in the (tan/3, Mjj± ) plane. 

5.2.2 Combined Fit 

We have performed a global Type-II 2HDM fit combining all the available observables (and using 
the tree- level SM predictions for the leptonic B decays). We find a global minimum Xmin ~ ^-^ 
Mjj± = 860 GeV and tan/3 = 7. Since the number of effective constraints varies strongly across 
the (tan/3, Mjj± ) plane, it is not straightforward to determine the proper number of degrees of 
freedom to be used in the calculation of the CL - even if the test statistic follows a distribution. 
According to the discussion in Footnote 7 on page 6 we avoid this problem by performing 2 000 
toy-MC experiments in each scan point to determine the associated p-value. The upper plot of 
Fig. 14 shows the 68%, 95% and 99% CL excluded regions obtained from the toy-MC analysis of 
the combined fit. For comparison the 95% CL contours using Prob(Ax^, njof) for n^of = 1 and 
"-dof = 2 are also shown. As expected, the n^of = 2 approximation is more accurate in regions where 
several observables contribute to the combined fit, while njof = 1 is better when a single constraint 
dominates over all the others (very small and very large values of tan/3) . For comparison the lower 
plot of Fig. 14 shows again the 95% CL excluded region obtained from the toy-MC analysis of the 
combined fit (hatched area) together with the corresponding regions obtained from the individual 
constraints. It can be seen that due to the increased number of effective degrees of freedom the 
combined fit does not necessarily lead to stronger constraints. 

The combination of the constraints excludes the high-tan/3, low-Mj:^± region spared by the B ^ tv 
constraint. We can thus exclude a charged-Higgs mass below 240 GeV independently of tan/? at 
95% confidence level. This limit increases towards larger tan/3, e.g., Mjj± < 780 GeV are excluded 
for tan/3 = 70 at 95% CL. 

''■^The stronger constraint obtained from the global CKM fit for B{B tv) is a result of the increased theoretical 
precision and, more importantly, the l.Qcr deviation with respect to the "tree-level" determination (c/. Table 5). 



5.2 Results and Discussion 



46 




Figure 14: Exclusion regions in the (tan/?, Af^±) plane. The top plot displays the 68%, 95% and 99% CL 
excluded regions obtained from the combined fit using toy MC experiments. For comparison the 95% CL 
contours using Prob(Ax^, ?T.dof ) for ridof = 1 and ridof = 2 are also shown (see discussion in text). The 
bottom plot shows the 95% CL excluded regions from the individual constraints given in Table 5, and the 
toy-MC-based result from the combined fit overlaid. 
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5.2.3 Perspectives 

Improvements on the low-energy 5-meson observables are expected from the KEKB and Belle 
upgrade program with an initial (final) target of 10 ab~^ (50 ab~^) integrated luminosity [132-134]. 
Parallel developments envision the construction of a new SuperB accelerator with similar target 
luminosities [135] . With respect to the 2HDM analysis, these programs are particularly interesting 
for the decays B ^ tv, B ^ fiv and B Dtv whose present branching fraction measurements 
are statistically dominated. Further improvement can also be expected for the measurement of 
B{B Xg^) with however less prominent effect on the 2HDM parameter constraints due to the 
size of the theoretical uncertainties. The measurement of the ratio of partial Z widths, R'^, could 
be improved at an ILC running at the Z resonance (GigaZ, cf. Section 4.3). The authors of 
Ref. [101] estimate a factor of five increase over the current precision, mostly by virtue of the 
increased statistical yield, and the excellent impact parameter resolution suppressing background 
from charm quarks. 

The LHC experiments will attempt to directly detect signals from charged-Higgs production, either 
via t bH^ decays, if Mfj± < mt, and/or via gluon-gluon and gluon-bottom fusion to t{b)H^, 
and the subsequent decay — > tu, or, if Mfj± > mt, via — > tb. The full tan/3 parameter 
space is expected to be covered for lighter than top (a scenario already strongly disfavoured 
by the current indirect constraints, especially the one from B{B — > Xg^)), while the discovery of 
a heavy requires a large tan/3, which rapidly increases with rising Mfj± [19, 20, 136]. 

6 Conclusions and Perspectives 

The wealth of available precision data at the electroweak scale requires consistent phenomenological 
interpretation via an overall (global) fit of the Standard Model and beyond. Such fits, mainly 
determining the top-quark mass, the Higgs-boson mass, the strong coupling constant, and the 
overall consistency of the model, have been performed by several groups in the past. The fit has 
sensitivity to confirm electroweak unification and the Brout-Englert-Higgs mechanism [137, 138] 
of spontaneous electroweak symmetry breaking for the dynamical generation of the fermion and 
boson masses, while posing problems for alternatives such as Technicolour in its simplest form [139], 
requiring more involved scenarios. Other theories, like Supersymmetry, are decoupling from the 
Standard Model if their masses are large. For such models the high energy precision data as well 
as constraints obtained from rare decays, flavour mixing and CP- violating asymmetries in the B 
and i^-meson sectors, the anomalous magnetic moment of the muon, and electric dipole moments 
of electron and neutron, exclude a significant part of the parameter space. However, the models 
can be adjusted to become consistent with the experimental data as long as these data agree with 
the Standard Model predictions. 

In this paper, we have revisited the global electroweak fit, and a simple extension of the Higgs sector 
to two doublets, using the new generic fitting toolkit Gfitter and its corresponding electroweak 
and 2HDM libraries. We have included the constraints from direct Higgs searches by the LEP 
and Tevatron experiments in the former fit. Emphasis has been put on a consistent treatment of 
theoretical uncertainties, using no assumptions other than their respective ranges, and a thorough 
frequentist statistical analysis and interpretation of the fit results. 
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Gfitter is an entirely new fitting framework dedicated to model testing in high-energy physics. 
It features transparent interfaces to model parameters and measurements, theory libraries, and 
fitter implementations. Parameter caching significantly increases the execution speed of the fits. 
All results can be statistically interpreted with toy Monte Carlo methods, treating consistently 
correlations and rescaling due to parameter dependencies. 

For the complete fit, including the results from direct Higgs searches, we find for the mass of 
the Higgs boson the 2a and Sa intervals [114, 145] GeV and [[113, 168] and [180, 225]] GeV, re- 
spectively. The corresponding results without the direct Higgs searches in the standard fit are 
[39, 155] GeV and [26, 209] GeV. Theoretical errors considered in the fit parametrise uncertainties 
in the perturbative predictions of Mw and sin^^^g, and the renormalisation scheme ambiguity. 
They contribute with approximately 8 GeV to the total fit error obtained for Mh for the stan- 
dard fit. In a fit excluding the measurement of the top quark mass (but including the direct 
Higgs searches) we find mt = 178.2 1|'2 GeV, in fair agreement with the experimental world av- 
erage. Finally, the strong coupling constant to 3NL0 order at the Z-mass scale is found to be 
as{M^) = 0.1193 lo QQ27, with negligible theoretical error (0.0001) due to the good convergence of 
the perturbative series at that scale. 

We have probed the goodness of the Standard Model fit to describe the available data with toy 
Monte Carlo simulation. For the fit including the direct Higgs searches it results in a p-value of 
0.22 lb 0.01_o.025 where the first error accounts for the limited Monte Carlo statistics, and the 
second for the impact of theoretical uncertainties (without these, the p-value is reduced by 0.04). 
The p-value for the fit without direct Higgs searches is similar (the reduced number of degrees of 
freedom approximately countervails the better value). The compatibility of the most sensitive 
measurements determining Mh has been estimated by evaluating the probability for a consistent 
set of measurements to find a single measurement that increases the overall of the global fit by 
as much as is observed in data, when adding the least compatible measurement (here ^pg). An 
analysis with toy MC experiments finds that this occurs in (1.4 it 0.1)% of the cases. 

We have analysed the perspectives of the electroweak fit considering three future experimental 
scenarios, namely the LHC and an international linear collider (ILC) with and without high lumi- 
nosity running at lower energies (GigaZ), all after years of data taking and assuming a good control 
over systematic effects. For a 120 GeV Higgs boson, the improved M\y and mt measurements ex- 
pected from the LHC would reduce the error on the Mh prediction by up to 20% with respect to 
the present result. The ILC could further reduce the error by about 25% over the LHC, and - 
if the hadronic contribution to a(Af^) can be determined with better precision (requiring better 
hadronic cross section measurements at low and intermediate energies) - a 30% improvement is 
possible. The largest impact on the fit accuracy can be expected from an ILC with GigaZ option. 
Together with an improved a(M^), the present fit error on Mh could be reduced by more than 
a factor of two. We point out however that, in order to fully exploit the experimental potential, 
in particular the anticipated improvements in the accuracy of M\y, theoretical developments are 
mandatory. If the Higgs is discovered, the improved electroweak fit will serve as a sensitive test 
for the Standard Model and its extensions. 

By extending the Standard Model Higgs sector to two scalar doublets (2HDM of Type-II) , we have 

studied the experimental constraints on the charged-Higgs mass Mh± and on tan/3, using as input 
branching fractions of the rare B decays B — > ^^7, B tv, B — > /iv, and B Dev, the Kaon 
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decay K — > /xz/, and the electroweak precision observable R^. Exclusion confidence levels have 
been derived by carrying out toy experiments for every point on a fine grid of the {Mjj± , tan/?) 
parameter space. At 95% confidence level we exclude charged Higgs masses Mjj± < 240 GeV for 
any value of tan/3, and Mjj± < 780 GeV for tan/3 = 70. 

Inputs and numerical and graphical outputs of the Gfitter Standard Model and 2HDM analyses 
are available on the Gfitter web site: http://cern.ch/gfitter. They will be kept in line with the 
experimental and theoretical progress. Apart from these update commitments, new theoretical 
libraries such as the minimal Supersymmetric extension of the Standard Model will be included 
and analysed. 
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A Standard Model Formulae 

This section gives the relevant formulae for the calculation of the electroweak observables used in 
the global electroweak fit. We discuss the scale evolution of the QED and QCD couplings and 
quark masses, and give expressions for the electroweak form factors and radiator functions. 



A.1 Running QED Coupling 

The electroweak fit requires the knowledge of the electromagnetic couping strength at the Z-mass 
scale to an accuracy of 1% or better. The evolution of a{s) versus the mass scale-squared s is 
conventionally parametrised by 

1 — /\a[sj 

following from an all-orders resummation of vacuum polarisation diagrams, sole contributors to 
the running a. Here a = a(0) = 1/137.035 999 679(94) is the fine structure constant in the long- 
wavelength Thomson limit [140], and the term Aa{s) controls the evolution. It is conveniently 
decomposed into leptonic and hadronic contributions 

Aa{s) = Aaiep(s) + Aa[Jd(s) + Aatop(s), (45) 

where the hadronic term has been further separated into contributions from the five light quarks 
(with respect to Mz) and the top quark. The leptonic term in (45) is known up to three loops in 
the q"^ ^ mj limit [141]. The dominant one-loop term at the Z-mass scale reads 



A„|-.,M|,.„ E (-^^I..f -.|.e,(^)) .3M...:o- ,4.) 

45 



^ Wz 

Adding the sub-leading loops gives a total of Aaiep(s) = 314.97-10"^, with negligible uncertainty. 



The hadronic contribution for quarks with masses smaller than Mz cannot be obtained from per- 
turbative QCD alone because of the low energy scale involved. Its computation relies on analyticity 
and unitarity to express the photon vacuum polarisation function as a dispersion integral involving 
the total cross section for e+e" annihilation to hadrons at all time-like energies above the two- 
pion threshold. In energy regions where perturbative QCD fails to locally predict the inclusive 
hadronic cross section, experimental data is used. The accuracy of the calculations has therefore 
followed the progress in the quality of the corresponding data. Recent calculations improved the 
precision by extending the use of perturbative QCD to energy regions of relatively low scales, 
benefiting from global quark-hadron duality. For the fits in this paper we use the most recent 
value, Aa[5^(M|) = (276.8 ± 2.2) • 10-^ from Ref. [66]. The error is dominated by systematic 
uncertainties in the experimental data used to calculate the dispersion integral. A small part 
of the error, 0.14 • 10~^, is introduced by the uncertainty in as{s) (the authors of [66] used the 
value as(M|) = 0.1176 it 0.0020 [142]). We include this dependence in the fits via the parameter 
rescaling mechanism implemented in Gfitter (c/. Section 3). 



''^While the two-loop leptonic contribution of 0.78 -10 ^ is significant (roughly one third of the uncertainty in the 
hadronic contribution), the third order term, 0.01 • 10~*, is very small, 



A. 2 QCD Renormalisation 



51 



The small top-quark contribution at M| up to second order in reads [143-146] to 

Aatop(M|) = -±^Mz fi^^_QQ2af 28.220 + 9.702 \n^)(a^\^i^))^ (47) 
45 TT mf I \ mf J \ J 



Ml 
mf 



0.1071 + 0.8315 af^ (/x^) + ( 6.924 + 1.594 In ] {^f (a*^)'^ 



0.7-10" 



(5) 

where the short-hand notation as = as /tt is used, and where as is the strong coupling constant 
for five active quark flavors, and jj, is an arbitrary renormalisation scale, chosen to he ji = Mz in 
the fit. 

The uncertainty on a(M^) is dominated by the hadronic contribution Aa^^^^M^), which is a 
floating parameter of the fit constrained to its phenomenological value. The errors due to uncer- 
tainties in Mz, mt and as are properly propagated throughout the fit. Other uncertainties are 
neglected. 



A.2 QCD Renormalisation 



Like in QED, the subtraction of logarithmic divergences in QCD is equivalent to renormalising the 
coupling strength {as = gg/An), the quark masses (m^), etc., and the fields in the bare (superscript 
B) Lagrangian such as af = s^Zagas, m^ = s'^Z^mq, etc. Here s is the renormalisation scale- 
squared, e the dimensional regularisation parameter, and Z denotes a series of renormalisation 
constants obtained from the generating functional of the bare Green's function. Renormalisation 
at scale fi introduces a differential renormalisation group equation (RGE) for each renormalised 
quantity, governing its running. All formulae given below are for the modified minimal subtraction 
renormalisation scheme (MS) [147, 148]. 



A.2.1 The Running Strong Coupling 

The RGE for asifi^) reads 



-Poal - Pial - I32a\ - f3^a\ - ... , 



(48) 
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The perturbative expansion of the /3-function is known up to four loops [149, 150] (and references 
therein), with the coefficients 
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(49) 
(50) 
(51) 

(52) 



where is the number of active quark flavours with masses smaller than ^u, and where ^3 — 
1.2020569. Solving Eq. (48) for introduces a constant of integration, A^"^^ with dimension of 
energy. The solution in the ultraviolet limit reads [151, 152] 
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where L = 21n(^/A("/)) > 1. 
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(53) 



As Us evolves it passes across quark- flavour thresholds. Matching conditions at these thresholds 
connect a^s^^ of the full theory with nj flavors to the effective strong coupling constant a^^ ^\ 
where the heaviest quark decouples. The coupling constant of the full theory is developed in 
a power series of the coupling constant of the effective theory with coefficients that depend on 
X = 21n(/i/mg) [151, 153-155]: 



l + Ci{x) [a 
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(54) 



with Os = Os/vr (recalled), and 
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(55) 



The integration coefficients Cj^o computed in the MS scheme at the scale of the quark masses are 



11 

C2,0 = 



82043 . 575263 2633 



C3,0 



+ 



27648" 124416 31104 



(56) 



The solution of the RGE (48) at arbitrary scale requires as to be known at some reference scale, 
for which the Z pole is commonly chosen. Three evolution procedures are implemented in Gfitter, 
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which lead to insignificant differences in the result. The first uses numerical integration of the RGE 
with a fourth-order Runge-Kutta method. The second (the one chosen for this paper) determines 
A^^) at Mz by numerically evaluating the root of Eq. (53), and the values for A^"/^^) are obtained 
via the matching conditions. Both methods use as(M^) as floating parameter in the fit. In the 
third approach, A^^) is directly determined by the fit without explicit use of as(M^). 



A.2.2 Running Quark Masses 

The MS RGE for massive quarks is governed by the 7-function defined by 

- 7(as) = -7oas - jia^ - 7203 - 7303 - • • • • (57) 



niq In fi"^ 



Its perturbative expansion has been computed to four loops [156] (and references therein), which 
for the c and 6-quark flavours reads [156] 

mcifi^) = 771^4^/^^ [1 + 1.0141 as + 1.3892 a| + 1.0905 a|] , ^^^^ 
nibin'^) = mbOs^/^^ [1 + 1. 1755 as + 1.5007 a| + 0.1725 a|] . 

The scale dependence of mg(/i^) is given by the scale dependence of as = as(/i^). The renormal- 
isation group independent mass parameters rhq are determined from the measured quark masses 
at fixed scales (c/. Table 1). 



A.3 Electroweak Form Factors 



The calculation of the electroweak form factors for lepton or quark flavours /, and k^, absorbing 
the radiative corrections (c/. Eqs. 11-13), follows the ZFITTER procedure [6]. It includes two-loop 
electroweak corrections [6, 16, 45, 157-160]. We use the intermediate on-shell mass scheme [6], 
which lies between OMS-I and OMS-II. These latter two schemes are used to estimate the un- 
certainty arising from the renormalisation scheme ambiguity (see [59] for more information). The 
form factors in the intermediate scheme are given by 
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(59) 
(60) 



where the superscript (G) stands for the inclusion of all known terms, whereas [G] = G + asG 
includes the electroweak one-loop corrections together with all known orders in the strong coupling 
constant. These QCD corrections are taken from [161]. The parameter 6p^^^ contains all known 
corrections to the Veltman parameter, defined by the ratio of effective couplings of neutral to 
charged currents [162, 163]. The subscript "rem" stands for "remainder". The correction Af 
given by 



[G] 
rem IS 



(61) 
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with the one-loop remainder 
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(62) 



where N^'''^^^ 
by 



3(1) is the colour factor, = sin^Ow and = cos'^Ow, and where Ap^ is given 
^P^ = -7^[^wwiO)-^^w{M^)] ■ (63) 



The terms S^^^(O) and S^^y(M^) are the W boson self energies discussed below. 

For the purpose of illustration we give the formulae for the one-loop corrections of the electroweak 
form factors at the Z pole for vanishing external fermion masses [164]: 
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where 
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[1 - 6IQ/I4 + 12Q24] ViziMl) + 
[S^^(O) - Sf,(M|)] . 



4- 
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The term S^^(M^) is the Z boson self energy. The vertex functions in the chiral limit are given 
by [16] 



Vivis) = ---2Rv -{3 + 2Rv)H-Rv) + 2{l + Rvf[U2{l + Rv)-U2{l) 
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(66) 
(67) 



where Li2 is the dilogarithm function, and where Rv = Ry — i'lv^ ~ ^v'^V / s, and Ry = My/s. 
The Z~7 mixing function in Eq. (65) is given by 



nf^(M|) = 2Y,Nl\Qf\vfU-s-m),m)) , 
f 



(68) 



where Vf = 1 — AQjs'^, and where the index / runs over all fundamental fermions. The integrals 
in (67) and hiQ^; Mf , M|) in (68) are given in Appendices C and D of Ref. [165]. 

For the two-loop corrections to the electroweak form factors, 5piem and dniera in Eqs. (59) and 
(60), the interested reader is referred to the original literature [6, 16, 45, 157-160]. Because of 
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missing two-loop corrections to the form factors and occurring in Z ^ bb, an approximate 
expression is used, which includes the full one- loop correction and the known leading two- loop 
terms oc mf. Non-universal top contributions [6, 166, 167] must be taken into account in this 
channel due to a CKM factor close to one and the large mass difference of bottom and top quarks 

where the function T^'^\m^ / M'^) is given in [166]. Since the first term in Eq. (69) represents one- 
loop corrections, it must be subtracted from the universal form factors to avoid double-counting. 
Let p'f^ and k[ be these corrected form factors (c/. Refs. [6, 166, 167] for the correction procedure), 
the form factors beyond one-loop are obtained by pb = p'l, (1 -|- rj,)^ and Kb = K'f,(l -|- Tb)~^. 



(69) 



A.3.1 Self-Energies of W and Z Boson 



The W and Z boson self-energies X^ww Szz ^^'^ on-shell derivative Yl'zz sums of 

bosonic and fermionic parts. The bosonic parts read [16, 168] 
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where the shorthand notation rw = Mfj/M^y and = Mfj/M^ has been used. The function 



L 



L{-s;M^^,M^^) is defined in Eq. (2.14) of Ref. [169]. 



A.4 Radiator Functions 



56 



The fermionic parts read [16, 168] 
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with rf = m?j/M'^ and Vf = 1 — AQfs\, (recalled from above), and where (...) in Eq. (74) stands 
for (— s; TTij^, m^^). The sums are taken over all fundamental up-type and down-type fermions of 
all SU{2) (g) U{1) doublets with masses mj^ and mj^, respectively. The integrals I„((5^; M^, M|) 
and J- are given in Appendix D of Ref. [165]. 
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The radiator functions Rv{s) and R\{s) absorb the final state QED and QCD corrections to the 
vector and axial-vector currents in hadronic Z decays. They also contain mixed QED Cg QCD 
corrections and finite quark-mass corrections expressed in terms of running masses. The following 
formulae are taken from [45]. They have been updated to take into account results from the recent 
3NL0 calculation of the massless QCD Adler function [17] (represented by the coefficient C04). 
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where the finite quark- mass corrections are retained for charm and bottom quarks only, i.e., all 
lighter quarks are taken to be massless. This restricts the validity of the above formula to energies 
well above the strange-pair and below the top-pair production thresholds, which is sufficient for 
our use. The mass rn'g denotes the other quark mass, i.e., it is frib if q = c and Trie l^ q = b. 
The running of the quark masses is computed in the MS scheme according to Eq. (58). The two 
parameters (5j{4) and 6co5 represent the next unknown coefficients in the perturbative expansion. 
They are treated as theoretical errors within the i?fit scheme, and vary within the bounds obtained 
when assuming a geometric growth of the perturbative coefficients with the perturbative order, 
i.e., for a coefficient H one has 6h„ = {Hn~i/ Hn^2) • Hn-i- 

The expressions for the fixed-order perturbative coefficients C^^'^'' in Eqs. (77) and (78) are given 
below. 

Massless non-singlet corrections [17, 170-173]: 
Co2^^-llC(3) + (-ii + ^C(3))n„ 
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which for Uf = 5 take the values C02 = 1.40923, C03 = —12.7671 and C04 
satisfactory convergence given that as(M|)/7r ~ 0.04. 
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Quadratic massive corrections [174]: 
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Quartic massive corrections [174]: 
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Power suppressed top-mass correction [174]: 
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Singlet axial-vector corrections [174]: 
37 7 

j(2) (^x) = - — + \nx + —x + 0.0132x2 , (100) 
12 81 

= -l^ + fc(2) + C(3) + |lnx + §ln2x. (101) 
Singlet vector correction [174]: 

R%) = (E^/) (-0.41317) a|(.). (102) 
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